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SYNOPSIS 


Thesis entitled, "Thermal Properties of Nearly and Weakly 
Ferromagnetic Fermi Systems: Spin Fluctuation Effects", 
submitted by SURESi-i G. kilSHRA 

in partial fulfilment of the requirement of the Ph.D. degree 

to the Department of Physics, 

Indian Institute of Technology, Kanpur 

J an uary 19 77 

The present work describes a systematic study of the 

effect of spin fluctuations on the thermal properties of 

( i) nearly ferromagnetic Fermi systems, e. g. Liquid Helium 3 

and ( ii) weak itinerant f erromagnets, i.e. ferromagnets with 

a low Curie temperature (T <<T F ) . Since the real (or virtual) 

c r 

transition temperature in these system? is much less than Tp 
(the basic characteristic temperature), one is always close 
to a phase transition for T<<Tp , and fluctuation effects 
are expected to be important. 

The temperature dependence of the spin susceptibility 
X and of the specific heat of a nearly ferromagnetic Fermi 
system (T <0) is usually discussed in terms of microscopic 
Fermi liquid theory or the paramagnon model. Results obtained 
in the latter calculations are valid over a very small temp- 
erature range T < Tp a , where Tp is the degeneracy temperature 

- 1 

of the free Fermi gas, and a is the Stoner enhancement 

o 

factor. One is thus unable to describe a number of significant 
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features of thermal properties for the degenerate range 
T<Tp. In the second and third chapters of this thesis 
we discuss a functional integral sch./in where the Hamilto- 
nian of a fermion system interacting via a short range two 
body interaction is transformed into that involving n 
interacting spin fluctuations ( n >_ 2 ) . This transformation 
is specially useful here because the 'Past’ electronic 
degrees of freedom are integrated out ? and the ’slow’ spin 
fluctuation degrees of freedom are made explicit. We show 
that the thermal properties of this interacting spin 
fluctuation system can be systematically analyzed. 

We analyze (in Chapter II) the contribution from 

one, two and three internal thermal spin fluctuations to 

the self energy of, say, the transverse spin fluctuation 

propagator (related *to X). In the low temperature regime 

(T < a Q Tp) , the paramaanon theory result is reproduced; 

while for temperatures a Q Tp < T < Tp a classical Curie like 

result is obtained. This latter behaviour is experimentally 

3 

observed in Liq . He^ in the corresponding temperature range. 

It is also shown that the spin fluctuation expansion is 

convergent, i.e. one need not consider the contribution of 

more than two internal thermal spin fluctuations to self 

energy. Results are compared with experimental results on 
3 

Liq. He , good agreement being obtained with experiment 
for T < 1°K (T/Tp'vQ.2). (This work has been briefly, re ported 
in Proc. 14th Int. Conf. on Low Temp. Rhys. ( LT 14) Vol. 1, 
p . 57. ) 
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In the third chapter C (T) is discussed. Retaining 
up to the harmonic (quadratic) term of the functional 
Hamiltonian, one obtains the RPA for the thermodynamic 
potential (and hence the specific heat). If the quartic 
term is considered in the quasi-harmonic approximation then 
one can relate the thermodynamic potential to spin fluctua- 
tion propagator. Using a quasi-harmonic approximation and 
the experimental results for spin susceptibility, ^ v (T)gp 
can be calculated. We have shown that noth the Boson and 
Fermion contribution (c.f. E. Riedel, Z. Phys. 210 , 403 

(1968)) are included in our approach. We get, for Q^(T)^p, 
a peaked curve with a linear rise in temperature in the low 
temperature side and with a slow fall on the h ig h ' temper a- 
ture side. In liquid He , for high temperatures (a Q Tp < T<Tp) 

a linear specific heat vs temperature relation is observed 

2 

experimentally. The slope of C /R vs x(=T/Tp) curve is ir /4 

2 

while for a free Fermi gas the slope is tt /2. This is 
accounted for by considering the suppression of density 
fluctuation excitations. This follows from the high sound 
velocity observed in Liq. He . A microscopic theory for 
this suppression is presented. A comparison of the theore- 
tical results with experiment shows good agreement up to 
TJ 0.2. 

In the fourth chapter, using the above calculations, 
we have discussed a T c = 0 ferromagnet. A classical Cyrie 



xi 

like behaviour down to 0°K is obtained. The lowest order 
spin fluctuation corrections are renorm aliz able and yield 
classical exponents. 

The next two chapters are devor.ed to the magnetic 
phase of a weak ferromagnet (T «Tp). The properties dis- 
cussed are magnetisation M(T) and spin wave stiffness B(T). 
These are usually studied in a molecular field theory 
( S to ner-Wohlf arth— Edwards ) or in RPA ( I zuyama-Kubo ) . These 
approximations seem to be quite inadequate in understanding 
the ferromagnetic problem as a whole, particularly at finite 
temperatures. A better approach has bean proposed by 
M ur at a-Do niac h , Moriya-K awabat a and R amak rishnan , who 
discuss the effect of mode-mode coupling on thermal proper- 
ties. In the fifth chapter, we have derived an expansion 
of free energy F in terms of the static, uniform z component 

E Z of the spin fluctuation field (< £ Z >^1'-'.). Spin fluctua- 
oo c 

~2L 2 

tion effect is included in the £ OQ term. The equation for 

the order parameter (obtained by minimizing ^(? o0 )) w hen 

simplified gives an Ed wards-Wo hi f art h (Jour. App * Phys. 

39 . 1 061 (1966)) Like result for low temperatures (T«T c ) 

2 

with a much larger coefficient for the T term. Near T c , 

2 

M , thus obtained;, decreases linearly with temperature. 

Spin fluctuation effects on Arrott plots are also discussed. 
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The penultimate chapter deals with an attempt to 
calculate D(T) in a self consistent and spin conserving 
approximation which goes beyond RP A anu includes spin fluc- 
tuation effects. A spin conserving approximation for the 

-U. 

transverse spin susceptibility X and for the related 

vertex function T' - is obtained by using Ward identity. 

The spin wave stiffness for small wave vectors can be 

obtained directly from the limiting value of a vertex y 

which is related to r + ~. One draws a sot of diagrams for i 

and analyzes various thermal c ontrib u cions. Roughly D«NI. 

2 3 /2 

We find for low temperatures T and T ** terms for D(T). 

2 

The coefficient oP T term is qualitatively different from 

that obtained by Izuyama and Kubo. The spin fluctuation 

5 /2 

interaction or 'magnon drag' term gives T dependence, 
as usual. We also show how the exchange splitting varies 
with temperature. (This work has been briefly presented 
in Mag. Letters _1_, 1 7 (1*9 76).) 


The final chapter contains a few concluding remarks. 



CHAPTER I 


INTRODUC TION 


3orutical study of ferromagnetism in a system 

of itinerant ferniigns' started with the work of Liloch, 

7 3 

51atcr“ and Stoner. Thoir work uses the Hartree Fock or 

molecular field approximation, where the two body electron- 

electron interaction is replaced by a self -co ns istently 

3 

determined one body potential. In the work of Stoner, 
the effective two body interaction is assumed to be of 
short range (on account of screening of the Coulomb poten- 
tial by conduction electrons). The temperature dependent 
properties of itinerant electron f erromagnets, such as 
the temperature dependent spin susceptibility x(T) in the 
paramagnetic phase, the Curio temperature T , the sponta- 
neous magnetisation H(T) have all been discussed in the 

4 

Hartree Fock approximation by Stoner, Wohlfarth and 

5 

others. The temperature dependence of physical quanti- 
ties arises in this model from the temperature dependence 
of the Hartree field, i.e. from the T dependence of 
fermion occupation numbers (Farmi Dirac distribution). 



2 


The dyn. r.,ic spin susceptibility ^(q,u>) and thence the spin 

, 6 
wave Dispersion h ve been considered in the RPA ( a level 

of approximation for response functions corresponding to 

Hartrec Fock appr oxim at ion (H F r,) for single particle 

energies) . 

The main thrust of recent work in itinerant fermion 
ferromagnetism has been in more realistic treatment of the 
ground state ( T =0 state), in particular, the effect of 
particle particle correlations , ^ The main qualitative 
conclusion of the work of Kanamori, Hubbard and Gutzwiller 
in this area is that the ’bare 1 short range repulsion U is 
replaced by an effective repulsion U_.p.p < U. This leads to 
more stringent conditions for the onset of ferromagnetism. 

The zero temperature spin wave stiffness has been calculated 

8 

by Callaway and Young. The temperature dependent proper- 

g 

ties can also bo calculated in this approximation, but one 
expects results qualitatively similar to that obtained in 
HFA, 'with U ff replacing U. 

In the mid-sixties, the role of incoherent spin 

fluctuations in nearly ferromagnetic fermi systems such as 

g 10 

Liq Hc^ and Pd was pointed out by Berk and 5chrieffer, 

11 12 

Doniach and Engelsberg and others. These particle 
hole excitations, known as paramagnons, are the low lying 
excitations of the systems, and affect strongly its 



3 


temperature dependent properties. In such systems the 

static zero temperature spin susceptibility x(T=0) is very 

large (:i.e. x(T=0)/Xp au j_i - a P >> 1 ) . It can be shown 

that the characteristic paramagnan energy is not Tp (the 

free gas fermi temperature) but a « T° Thus, as 

r o r 

shown by leal Nonod et al . „ 1 

X(T) = X ( T =0 ) - A( T/Tp a Q ) 2 ; 

i.e. x(T) varies rather strongly with temperature, on a 
scale T° a . Because of this, the general impression has 
been that paramagnon theory results are valid in the rather 
narrow temperature range T<<Tp a . 

A somewhat different aspect of spin fluctuation effects 

14 

was first discussed by Doniach and hurata. Using a func- 
tional integral transformation, they showed that the Stoner 

Hamiltonian (H = ) £. n + U j n, (r) n.(r) dr) could be 

k,5 k ko * + 

transformed approximately into a Hamiltonian for classical 

interacting spin fluctuations m . Assuming that the number 

of classical spin fluctuations (i.e. those with energy <<kgT) 

is comparable to the number of degrees of freedom of the 

system, they obtained in a Hartree like approximation for 

fluctuations, a Curie-Weiss law for x(T). A very different 

1 5 

analysis by Mariya and Kawabata is based on a self- 

consistency requirement on x(„Q) • They assume that the 

-1 

RPA expression for x(.s) is modified by an additive 
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constant. This constant is determined by requiring that 

2 2 

the susceptibility satisfies the identity x = (9 F^/SFl ) jV ,^_j, 
the free energy being determined as the sum of ladder and 
bubble diagrams. This work st abl ished that near T , the 
number of thermal spin fluctuations is equal to i\l (number 
of degress of freedom of the system) x (T/T°)^. Moriya, 
Kawabata and coworkers 10 have since then discussed spin 
fluctuation effects in magnetisation, spin lattice relaxa- 
tion time, in itinerant electron antif erromagne ts . There 

1 7 

is experimental evidence for the effects calculated. 

1 8 

Ramakrishnan attempted to discuss spin fluctuation 
effects from first principles, starting from an equation 
of motion for x(g> w ) which leads to a three particle 
Green's function G*^, making a diagrammatic interacting 

spin fluctuation analysis of the latter. He showed how 
spin fluctuation effects could be decomposed into zero 
point and thermal parts. Further he found that the 
coupling constants, such as the critical coupling for 
ferromagnetism, and the four fluctuation coupling constant 
were likely to be very different from their zeroth order 
value because of the intermediate coupling nature of the 
problem. He further estimated fluctuation correlation 
effects which, very close to f , dominate the mean fluc- 
tuation field term (responsible for a Curie-Weiss law). 

More recent papers on the subject are due to Hertz and 
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'I u 2 G 21 22 23 

Klenin, '' nurata, K aw ab at a , " Yamada, Hertz, Gumbs 

and Griffin. ^ 

Inspito of the large amount of work in this area, 
there has been no systematic attempt to discuss spin 
fluctuation effects in nearly and weakly 

ferromagnetic ferrni systems from first principles. We try 
to do this in the present work. 

The area of our concern can be described by the 

parameter space of Fig. 1, where the x axis represents 

- 1 

the short range interaction U in the units of P £ p. The 

y axis is temperature k T in the same units. The fermion / 

B / 

system has a paramagnetic ground state for small U, and 

presumably (i.e. if it is not a hott-Hubbard insulator, 

ant if erromagnet ) a ferromagnetic ground state for large U . 

The boundary is located at ijp £ p ,= 1 in HFA. Earlier work 

on correlation effects can be viewed as an attempt to 

rescale the x axis, i.e. to find ( U , P £ ) and to Inc ate 

the point U ^ P £ p = 1. Starting from this point, a curve 

can be drawn in the (T,U) plane, separating the paramagnetic 

and ferromagnetic phases. Close to U P £ p = 1, the curve 

2 

(marked SW) is parabolic in the HFA, i.e. (T^ Pep) a (UPgp-1) 
This follows simply from the nature of the degenerate ferrni 
distribution. In the work of Koriya ct al . this is 
replaced by a nearly straight line marked ISDIHK. Consider 
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the shaded region |Up £ p-l| << 1 , T P £ p« 1 . The system is 
always degenerate, and since (for dp £ pl>, 1) m^ (saturation 
magnetization at T =0 , in dimensionless units) a ( U Pgp -1 ) <<1 , 
magnetic ordering affects electronic energies only weakly. 
Clearly, in this region, one is always close to the ferro- 
magnetic instability. Thus longitudinal and transverse 
spin fluctuation excitations ai ways have characteristic 
energies much lower than Pep, so that the temperature 
dependent properties are expected to be dominated by these 
excitations. 

The nearly ferromagnetic regime is the region A, 
bounded by the line Up £ p=1. in the area close to kgTp £ ^ = Q j 
i.e. in the lower half (bounded by Tp £ p =: ( 1 -U p £ p ) ) lies 

the 'paramagnon' regime. In the part A^ one is probably 

in a classical spin fluctuation region, as is also plausible 
from its contiguity with 3, the paramagnetic regime, where 
classical spin fluctuations dominate the temperature 
dependent properties. The weakly ferromagnetic region L 
is also divisible broadly into a 'classical' regime 
close to the transition line and a region analogous (in 

its characteristic temperature dependence) to the 'paramagnon' 
regime A 2 . 

In this thesis, we first consider the region A. We 
show how the Hamiltonian of the system can be transformed 
into that involving n interacting spin fluctuations ( n > 3 ) * 
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The transformation maintains rotational invariance and no 
'classical' or zero frequency approximation is made. The 
fast electronic degrees of freedom are integrated out, and 

l 

determine the s.f. coupling. Uc first consioer the X( T ) of 

this interacting s.f. system, and analyze the contribution 

from one two and three internal s.f.s to the self energy £ 

of say the transverse s.f. propagator (related directly 

to x)* Such a spin fluctuation interaction expansion is 

shown to be convergent; i.e. the leading terms in £ involve o:ss ^ 

and two correlated internal s.f.s. The three s.f. contribution 

is less by at least one power of T p £ p (see Chapter II for 

details). Further, it is seen that the same physical 

process (diagram) is called the paramagnon process in the 

low temperature or regime and classical s.f. in A^ . 

One thus obtains in a single scheme, a theoretical express iaT\ 

for x(T) which gives the leading term for the entire region 

25 

A. We compare our prediction with experimental results 

2 

for X(T) of liquid He . He is a pure homogeneous fermi 
liquid, and thus is free of unknown complicating effects 
due to band structure (density of state structure) present 
in nearly magnetic metals. Further, by increasing pressure, ■■ 

one can come quite close to the boundary UPep=1. We find ; 

very good agreement with experiment, over a temperature > 

range 0.0 < t = T p £ p < 0.2. . 
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We then study (in Chapter III) the thermal proper- 
ties of s nearly ferromagnetic fermi system, and calculate 
the specific heat in an approximation which is shown to 
retain the leading term. Again comparison is made with 
liquid He , for which C (T) data are available. Good 
agreement with experiment is obtained over a temperature 
range 0 <x<0.2. The shape of the C V (T) curve showing a 

low temperature bump and a linear flattening out C V (T) - 
TT 2 T 

— -- CTT ( i . a . half the free fermi gas value) are both explained. 
4 p 

The former is due to spin fluctuations and the latter to 
the facx that density fluctuation excitations have a high 
characteristic energy (compressibility 'vl/IO of free fermi 
gas compressibility). 

We believe that the very good agreement between 
theory end experiment over a broad temperature region argues 
strongly for the quantitative correctness of our systematic 
theory, since He is a very clean case. bupporting evidence 
is present in the experimental X(T) values of various nearly 
ferromagnetic metals and alloys, e.g. iMi^Ga, HfZn 2 , NiRh, Hd, 
YCo 2 , U 2 C 3 etc. Their x( T ) curve is quite similar to that 

3 

discussed above for He . 


The region B, i.e. the paramagnetic region, has been 
investigated by several authors, as mentioned above. Here, 
the one and two internal s.f. approximation for S becomes 
insufficient close to T . An estimate of the three s.f. 



contribution, made by R am akris hn an , can be used to roughly 

demarcate the critical region (Fig. 1, ohoded area close to 

the line hDMK). An interesting point is the shrinkage of 

the critical region as T c -»• 0 . In particular, a T c =0 ferro- 

magnet is expected to have purely classical (mean field) 

critical indices. This has also been painted out recently 
23 

by Hertz using a renormalization group method. In 
Chapter IV, we present explicit calculation for X( T ) , 
confirming this expectation. 


The remainder of the thesis is devoted to the 
region C, i.e. the fe rroma gne tic all y ordered phase of a 


weak ferromagnet (m <<1 or T p <<1). We obtain an 

O C b. p 

2 

equation for m(T) in an approximation scheme retaining 

the leading thermal spin fluctuation effect. We find 

that m(T)^ is strongly coupled to spin fluctuations. 

Longitudinal spin fluctuations have low characteristic 
2 

energies m £p. The transverse spin fluctuation spectrum 

2 

has low lying spin wave pole at frequencies oj =Dq ( d “ m, 

* q 

roughly) and has a resonance structure in the Stoner conti- 
nuum. We calculate the spectrum of those fluctuations and 

2 

find their effect on m . The effect is seen to determine 
its temperature dependence; that arising from temperature 

dependence of the fermi distribution is very weak. At low 

2 2 2 ( T p ep ) 2 

temperatures (T p « m ), we find m (T)=r.i ( 1 -A j ) 

C'p q m H 

‘ o 


( 1 -A 



where A is a constant. The spin wave contribution to 
2 

m (T) is smaller by a factor rn . This is in accord with 

o 

2 

experiment. Well above T /2 , the T dependence changes to 

2 

a classical form m <x (T -T). We do not attempt a detailed 

comparison between experiment and theory for obvious 

reasons. However, the study of magnetisation M(H,T) in a 

magnetic field H provides a case where a qualitative new 

2 6 

prediction can be made. Here, Wohlfarth et al . have 
shown, by analyzing a Ginzburg Landau like model with a 
5 toner-Wo hlf arth temperature dependence for N(T), that the 
Arrott plot of M^(H,T) against H/h(H,T) should be a straight 
line. The same paper exhibits Arrott plots for IMiPt alloys; 
these are systematically and decidedly curved for small H/l v l. 
We obtain the Arrott plot expression in an approximation 
scheme retaining spin fluctuation effects. We show that 
the deviation from linearity can be qualitatively 

explained for small H/M and large temperatures. 

In Chapter Vi, we take up the question of the tempera- 
ture dependence of the spin wave stiffness D in a weak 
ferromaynet. In the RPA,^ D(T)“ m(T) for small m; this 
m being the spontaneous magnetisation in the HFA. Beyond 

the RPA, there are very few results. A general phenomeno- 

2 7 

logical analysis due to Izuyama and Kubo leads to the 
low temperature form 
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D ( T ) = D ( T =0 ) i AT 5/2 - BT 2 , 

the coefficients A and B being unknown. This is also con- 
firmed by the microscopic analysis of de Pasquale and 
2 8 

Corrias for a strong ferromagnet ( = N, =0). We 
arc interested in knowing whether and how the spin fluctua- 
tion excitations affect D(T). For this purpose, we start 
from an approximation for electron self-energy X- which 
includes the effect of electron scattering from spin 
fluctuations. Then, by fallowing a standard method, we 
construct a spin conserving approximation for the transverse 
or spin flip response function. From this we determine 
the spin wave stiffness D. We find that at low temperatures 

D(T) = D (1 - A' T 3/2 „ B'T 2 ) , 
o 

3 / 2 

i.e. there is a T term due to transverse spin fluctua- 
tion (spin wave) contribution to single particle self- 

2 

energy, and a T term arising from longitudinal and 

incoherent transverse spin fluctuations. The size of the 

3 /2 2 

T term, however, is smaller than that of the T term 

by a factor m. There are no experiments on i)(T) in weak 

ferromagnet s , to compare with our prediction. However, 

2 

large T terms are reported for D(T) in many itinerant 

2 9 

electron ferro magnet s. 

In the concluding chapter, we point out a few areas 
where further work is needed, and also other phase transi- 
tions in fermi system where such fluctuation effects might 
be significant. 
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1 . Intro duction 

In this chapter we address ourselves to the problem 

of the temperature dependence of spin susceptibility x( T ) 

30 

of nearly ferromagnetic ferrni systems. These systems are 
characterised by an enhanced zero temperature suscepti- 
bility. An ideal example of such a system is normal liquid 

31 3 2 

helium-3. There are some metals, e.g. palladium and 

,, , r 33 u _ 3d ,, r 35 vr 36 .,. 0 .37 

alloys like iMi^Ua, HfZnjs, U 2 L 2 » YLo^* n in h 

which also show such behaviour. In the calcula- 

tion we will refer to liq. Ho only, as there are no band 

structure and density of states complications. We can, 

2 2 

therefore, work in the parabolic band model e^=1i k /2m 
and with free fenni gas density of state, p( £ ) . 

The most recent and precise experimental work on 
3 2 5 

X(T) of L s.q He is by J. R. Thompson and his associates 
Duke University. The work reveels two characteristic 
features in two different temperature regimes. 

( i) At very low temperatures the data can be 
represented by the expression 



1 4 


X(T) 


i. 

a 


(1 


-r O 

F 


•) 2 > 


a 


( 1 . 1 ) 


1 his douoly enhanced T dependence has been predicted by 

paramagnon theory, bur the oi. served behaviour is over a 
2 

wider ( x behaviour even for x 'v a ) temperature range. 


(ii) At sufficiently high temperatures (T > 0.5°K) 
the susceptibility is found to tend asymptotically towards 
the Curie law relation, with a large Curie constant. dote 
that though the system is highly degenerate (T << T°), its 
magnetic behaviour is like that of a collection of inde- 
pendent classical moments. This is surprising. 


A good microscopic theory should explain these 
features, at least qualitatively. We direct our effort 
to attain this goal. Before doing that, let us pause for 
a while and review earlier attempts in this direction. 


There are sevcr 0 l theories with which the 
experiments in temperature range T << O.S^K can be compared, 
For an ideal farmi gas with a degeneracy temperature 
T° = --- (3Tr n) 7 /k , the susceptibility )(_ can be 
represented, for t < 0.1 by the expression 

Xid 


Xp (1 - 12 T 


...) 


( 1 . 2 ) 


For LIq He^ the effect of interaction between spins has the ef 
feet of increasing *% Q over the expected value by a factor vary 



between 10 and 2G according to t hu density. The Stoner 


3 


theory for a fermi gas with interaction then predicts 


• 1 


-1 


U , u is the density dependent quantity prop or- 


x = x id 

tionai to the exchange integral. Assuming U to be tempera- 
ture independent, one obtain-., in the low ten, per ature limit 


*C 


X P 


= a. 


(1 - 


2 2 
IT T 

12 a o 


(1.3) 


a singly enhanced T dependence. The next atte rnp t w as 

3 8 

that dun to Richards who calculated lowest order tempera- 
ture corrections to fermi liquid theory result and expres- 

2 

sad the coefficient of T term as function of certain 
parameters obtained from a combination of specific heat 
and susceptibility results. 


The most remarkable result is due to Beal No nod 
et al . , who calculated X(T) near T=0 r K in the paramagnon 
model. They calculated the correction to free energy due 
to transverse and longitudinal paramagnons (or spin fluctua- 
tions) represented by a collection of ladder and ring 
diagrams respectively. The corr esponc ing correction to 
zero 1; ci.,p-.r at ure susceptibility is found for T << to be 

X (T) = — {1 - (- a T -) 2 + 3 {(x/a ) 4 In (x/a )}> (1.4) 

a 2 4 cl o o 

a o 

There is no adjustable parameter introduced and the magni- 


tude of the third term is small in comparison to the second 
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one for t<<a Q . Also notice that there is no T In T 

term in the susceptibility, Cb u ch a term is obtained by 

3 5 4,1 

usrnea and Mis aw a. ) The theory works quite well in the 
10W r a i.ip or at ure regime * It appears as though the theore- 
tical expansion parameter is (x/a ) whereas the form (2.1.4) 

o 

seems to represent experimental results accurately up to 

x /a * 0.3(7). This app ro ac h is also m arr ed by a 1 ac k of 

self consistency, first pointed out for itinerant electron 

1 5 

ferromagnetism by f-ioriya and Kawabats. 

In the region of 0.5 C K and above where the suscepti- 
bility tends towards the classical Curie like behaviour, 
the theoretical prediction of X(T) is a difficult problem 
because the thermal energy k n ,T is higher than spin fluctua- 
tion energy £p and lower than the fermi energy kgTp and 
the interaction energy U. The f iur at a-Do niac h approach 
involving classical interacting spin fluctuations, has been 
successful in predicting a Curie-Weiss law for X ( "O for 
magnetic systems (T > T ). 

However, sc far there exists no microscopic theory 
applicable over the entire temperature range, i.e„ one 
which describes accurately both the par amagnon behaviour 
and the ’classical' behaviour, as well as the intermediate 
region, with well defined limits of validity and accuracy. 

We consider, in this chapter, a microscopic spin fluctuation 



theory end show that it is prcsible to understand the 
temperature dependence of X~ 1 ( T ) over a wide temperature 
range t << 1 (but not necessarily x/a << 1 ) . The leadinq 

Q 

-1 3 

terms are then used to calculate x (T) for He , for T < 0,2. 
The agreement with experiment is good, and the deviation 
at high t explicable as to size and sign. We summarize our 
theory now, details being given in succeeding sections. 

Obviously, the low lying ( | q" | <<k~ ,co <<Ep) spin 
fluctuations ( paramagnons) are important in determining 
the low temperature (x<<a Q ) behaviour of the system. Their 
effect dominates over that due to uncorrelated or 'free' 
particle hole (Stoner) excitations because of the exchange 
enhancement. It is. suggested that if the effect of inter- 
action between spin fluctuations is calculated self 

consistently (mean fluctuation field approximation, MFFA), 
one may >,st a qualitative understanding of the behaviour 
in the high temperature (x>>a but x<<1) regime also. In 
the functional integration scheme, u.. t 0 xiny the trace 
over the fermion variables, we have an effective inter- 
acting spin fluctuation Hamiltonian F ( £ ), with n spin 

q s u 

fluctuation coupling (n >_ 3 ) . The 'bare' coupling cons- 
tants can be evaluated in terms of the free fermi gas 
parameters and depend mildly on momentum (q/2k^), energy 
(to/e.p). The fluctuation propagator 



is directly proportional to the spin susceptibility func- 
tion X ae ( <?, z )n ) of the system We therefore analyze D. 

The bore fluctuation propagator corresponds to RPA. We 

are interested in the temperature dependant contributions 
-1 - 1 

to D (jj - D_ ~ Z) . These can be isolated in principle 

G 

by performing the frequency sums in diagrams for E , and 
considering terms with one or more positive frequency 
bose p or ts (i.e. 0(w) (exp ($w) - 1)“ ). Introduction of 

one internal spin fluctuation line in Z leads to the 

, . 30 

2 qu at 10 n 

X(q) X°(.q) _1 = (1 - d eff X(q) - X ^XU’))" 1 U.6) 

for susceptibility. The re nor hi ali zed four fluctuation 

coupling vertex X T cannot be evaluated reliably. because 

of its mild dependence on co ? q anc * T, we approximate it 

by a constant Contribution from incoherent fluctuation 

terms and zero point (T=0) part of X term are lumped 

together and are assumed to renormalize U (to . For 

T << (X ✓ the equation (1.6) need not be solved self consis- 
o 

tentiy and one has the par ainagnon theory results 

* n 

A (X" ) « T^/a . F or temperatures a Q < r <<G "the self 

consistent solution loads to a Curie like behaviour; 

, -1 . 

(X ) “T . 

The above expression for X i s useful if correction 
terms are small. The next thermal correction terms for Z 
are due to two and three internal, correlated spin fluctuations 



Qualitatively, the contribution of Z^ (correlated) is 
similar to one s.f. term,. The resulting contribution to 
y ' Xn from 1 is estimated to be less than x^ os - ^ for 

l □ 

T < a end thus vary small. f jr t << I but T > a , the 

Q w — C 

n 

estimate is x~ In (1/a), T hr mean particle field term ^ocs 
1 T ~ . Thus these corrections arc of relative order x, 
and the solution of the above aquation is expected to give 
a correct X ( T ) fort << 1 . This smallness of correlation 
effects is due to only a small fraction of the total 
degrees of freedom being excised since T « 1°. 

Wa have made a detailed calculation of tho suscep- 

3 XI 

tibilifcy of Liq He , using equation (2.1.6) plus 2 

corrections also included. Thu agreement is very good 

up to ^ = 0.2 and becomes progressively poorer. The 

fluctuation correlation correction is of the correct sign 

and size to explain the difference between theory and 

e xp er imo nt „ 

To summarize, we have a scheme where the effect 
of one, two and higher thermal spin fluctuation can be 
calculated systematically from the first principles. We 
have shown that if the paramegnon effects are calculated 
self consistently, one can understand the temperature 
dependence of spin susceptibility over a wide range of 
temper ature . 
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In the fallowing sections we derive the functional 

Hamiltonian >. nd calculate the cPFect of one, two and three 

internal, correlated thermal spin fluctuations. Thun we 

compare our results with the previous theories and the 

3 

experimental results on Liq H: and other metals and alloys. 

The calculations and analysis presented in this chapter are 
rather detailed and will be referred to quite frequently- 
in the subsequent chapters. 

2 . Fluct uation Interaction T heory 

It is possible to transform the basic Hamiltonian 
of the system, given in terms of the Fermion variables, to 
an effective Hamiltonian expressed as a power series in 
interacting fluctuation fields (Boson variables). This 
is done by using the well known method of functional integra- 
tion. The re sult is that the slow spin fluctuation 

degrees of freedom are made explicit, their 'bare' spectrum 
and the coupling between them being determined by a trace 
over the fast electronic degrsos of freedom. Clearly the 
trans Formation is particularly useful if oj s ^ <<ai f aS t* 

In our case this translates into T c << Tp for the weakly 

- 1 - 1 

ferromagnetic case and x(T=0) <<'Xp a(J j_^ for the nearly 

ferromagnetic case. We analyze the effect of some lowest 
order terms in such expansions. 
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2 s . F unc tiona l Integral F or mal ism 

Let us start by discussing the functional integral 

Formulation of the theory of itinerant Fermion systems. The 

45 

rules tor such a transformation have been given by Rice. 

We simply transcribe them for our case. Here the statistical 
operator is represented as a functional integral (a Gaussian 
functional average) over a 1 tirnn 1 ordered exponential, 

Feynman time ordering is necessary when different terms 
in the Hamiltonian lin the exponent ial) do not commute with 
each other. 

The Hamiltonian of an interacting fermion system 
is given by 

l 5^.5+ (2a. 1 ) 

■$ 9 

For detailed discussion about this form the appendix I may 

be referred. Hero and ar* are annihilation and creation 

k k a 

operators For fermions of momentum k and spin cr o is 

the single particle energy measured from the chemical poten- 

*"■4 

tial V. o->- and S-* arc Fourier transforms of number density 
q q 

and spin density operators respectively. We have written 
the interaction in a form which is explicitly invariant under 
rotations. An effective zero range interaction between the 
fermions is parametrized by V and U which can be related 
to experimental properties through dynamical models. 
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Thu grand partition function is defined by 


t r a x p [ - 0( H- pi-0 3 


(2a. 2) 


where 0 = 1 /k, ; T . Substituting h- Pit* from the equation ( 2 a. 1 ) 
and introducing an ordering label s, we write 


2 = tr exp { - ^7 [ l 


E k "kos d kp 


a, + 


V 1 P 


qs -qs 


5 0 . s ] } 

os -os 


( 2 a. 3 ) 


The sum ovsr s runs from 1 to k. With the introduction of 
the ordering label, the operators in the exponent can be 
treated as c numbers with an error which goes to zero 
as id -»<*». Using the Str atonovich-Hubb ard identity, 


:p { |a| } = / 


oo dx 1 dx2 


exp {- I x I 2 + ax* + a*x} (2a. 4 ) 


where x = x ^ + ix 2 etc., one can express the partition 

function as a functional integral over boson variables f, £ 

df d | . 

A = L q . s , i =1 ,2 .A P 4 T J 

"*P <-Uu= 1.3 (if q S ' 2+ 'ti' 2 ’ * 


x L (f 


(2a. 5 ) 


where 


L t f qs . I qs i = tr e *p { - •; bl e k B 


k kps kps 


- y C p ( p f* e + h,c . ) 

£ v K qs qs 

qs ^ 

- Y c S (S j ^ _ * + h.c. )1 } 


. y c S (S J * + h.c 

L . qs c s 

qsj 


j — z , + , — . 


( 2 a. 6) 
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1 it h c^ = (-£) 


VI/2 


s - 


(2a* 7) 


Hera L behaves like the partition function of a system of 
non- in Lur acting fermions moving in presence of u set of 


random time dependent fields. As i\! 


the set of points 


-7- 

{5 S ) etc., becomes a continuous function £( u ) defined on the 

interval [0,3] with the boundary condition 5(0) = 5(8) etc. 

- f 8 

and 8/N L J du. The above condition originates 

s o 

from the boundary conditions imposed on Green's function 
and operators in the imaginary tine domain. Wow the 
quantities depending on u can be expanded in a Fourier 
series as is usually done in finite temperature many body 
theory, and Eqn. (2a. 5) becomes, 


exp ( - l ( | f a | 2 + I 5 Z 1 2 + 1 5^ I 2 + i 5“ ' ; 

a 

, 5” } (2 S .8) 


df 

a 

d E . 
ai 

11 • 

a> i 

V 


L tt, * 

5 Z 

, 5 + 

a 

a 


a = ( q , 

z ) 
m 

and 


L tf a • ?a > * T u e 


3 4 - 

xp {- / du Ka (u) a (u 
■ o k a 


t c<> f B du t (f* P + a 2,,i " uj/B + h.c.) 

8 ‘a a a W 

z 1 P , r , z* r Z 2fimu/3 l_ _ \ 

+- c I d u ) ( F b e + h.c . ) 


+ C z / du I ( ? Z b e 

8 l a qu 


, v / t- + 2 rnimu/ 3 , , ^ \ 

J du £ ( £ 3 c + h . c . ) 

o a qU 

f e du l (5~* S“ 0 2uim ^ + h.c.)} (2a. 9) 

J L qu 

o a ^ 
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to expand Y[f,£] in powers of the fields, i. 


^ ~y -y ~y +y 

Y [ f * 5 ] = +1 K 2 ( q, m) | f o m | + y K 4 ( q 1 , m <j ; q 2 , m 2 ; q 3 , m 3 ; , m 4 ) 

n m * m r 


q ? m 


q . 5 m . 

1 1 X 


f(q 2 m 2 ) f (q3v m 3) f(q 4 ,™ 4 ) 5 ( l q ± ) ( I m .j_) + •• 


i = 1 


a,{> ->* | r agy6,“> -*• , 

l b 2 ( ^ m) ISqJ + +1 ... b 4 (q 1 rp i ? q 2 m 2 ° q 3 m 3 * q 4 m 4 ) 


L 

qi.t 


C: . , m . 

1 1 X 


«( q | ) %( q-> nr,) ^ y( " q- 1.I3 ) ^<S(q 4 m 4 ) 5 ( I q i )5 ( \ rr^ ) 

i=1 i=1 

-i- . ... + terms involving coupling between spin and 

density fluctuations; (2a. 13) 

where a, 3, y, <$ are the spin indices. The coefficients K's and 
G's are given in terms of closed fermion loops. For example. 


r a $ 


(q,m) = (1 - U ^(q,m)) 


where y is the functionaJ. evaluated by Lindhard for a free 
''■a 

fermion model. In the interacting fluctuation fields, these 
coefficients can be regarded as the vertices. 

Clearly the form is analogous to that discused oy 
Ginzburg, Landau and Wilson for classical interacting 
fluctuation fields. The quantum effects can be traced to 
the non-commutativity of H q and H ' ., which force us to 
write cxp(-gH) in the interaction representation, requiring 
the functional averaging identity (2a. 4) to be applied 
for each time u. This makes the order parameter time 
dependent, with the consequence that (hatsubara) frequencies 
(z ) appear on the same footing as wave vectors. 
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bow we can write the interacting fluctuation Hamil- 
tonian m a physically motivates approximation scheme and 
carry out the functional averaging. This gives the partition 
function ana hence the thermal properties . In the next 


chapter wc 

will 

use this procedure 

to c ale ul ate 

the so 

ec if ic 

heat of L:.. 

H, 3 . 

Here wo consider 

the magnetic 

s u s c B ; J 

ti- 


bility arid its tcjmperature dependence. From the expression 
for { [f, it is seen that the thermal properties are 
determined by the density f luc tu at ions., the spin fluctuations 
and their interactions. Since the contribution of the 
interaction of density fluctuation and longitudinal 
spin fluctuations to X(T) is quite small, we neglect it. 


An expression for the static spin susceptibility 


which is defined as 


1 r In Z 

3 


(2a. 1 4) 


h * u 


can be written in ter;, is of the functional average. As the 

field n enters additively j.ith in j_, a simple change 

■? 2 gi ■ h 

of variables from to c, -• — yields the expression 

00 00 z J r 


of variables from E~ to c 

00 GO 

y ( ) ,2 

X (0,0) = JL [ 2 < 5 Z > 

U 00 


( 2 a. 1 5 ) 


where < > denotes the functional average. The above result 
can be generalized to that for the dynamical susceptibility 
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™ 0 -y. 

where D M (q, z^) is the fluctuation 'propagator'. Thus to 
calculate the response of the system we have to calculate 

9 2 

the averages <|£ Z f> and < | £ + { >, etc. 


2 _ 2 _ 2 _b_ J5,p_i n_ F luctuat io n Ex p ans ion 

The fluctuation propagator U c an be obtained from 
Y t?3 ? using the diagrammatic rules given in Appendix II. 
The bare propagator is, as usual, the inverse of the coef- 
ficient of the quadratic term in the Hamiltonian, 


if =11 - U X°(a)l (2b. 1 ) 

c 

Denoting the totally free fluctuation propagatox (Li=0) by 
a dotted line, , the above equation can be expressed as 



where the wavy line represents the propagator D q , and the 
continuous lines are fermion propagators. This approxima- 
tion, where no fluctuation interaction effects are included, 
is the riP A . b 


The lowest order self energy diagrams for the trans- 
verse s.f. propagator which result on including the quartic 
interaction effects era shown in F igs. 2a to 2c. The propa- 
gator D is given in terms of 2 by the Dyson equation 




o 


(2b. 3) 
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(2g) 


F ig. 2 . 


Self energy diagrams for the transverse spin 
fluctuation propagators 
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Highur order diagrams are sho-;,n in Figs. 2d to 2f # Wow we 
systematically analyse the t j mper ature dependence of the 
contribution due to these diagrams. 


The temperature dependence* of E(q = u . m = u) arises 
from internal frequency summations which can be connected to 
integrations over frequency with a Bose factoi, i.u. to 
/ <j>(oo,T) (oxp (|3io)-l)~ „ This can be divided into a zero 

— 03 

point port and a thcjrmal part. The former is non zero at 

zero temperature and depunds very weakly on temperature 

00 

through the explicit temperature dependence of J 4>(a) y T)dw # 

— CO 

This function depends explicitly upon the fermion variables 

only and can be shown to bo characteristically of the form 
o 2 

A + d ( T/Tp ) at 1 o w temperatures. Wo omit such terms. The 
thermal contribution is expected to bo sizeable if there arc 
considerable low energy fluctuations, and if their coupling 
is strong. The former is true near a ferromagnetic instabi- 
lity and the latter is true in our case because the coupling 
is via electrons. Thus the coupling nac t size U , range 'v 
2k.p and in wave vector an • energy space respectively. 

Consider first the thermal part of the diagram 
(Fig. 2a). Its contribution to Z is 


Z 2a ( 0,0) T"g G 1° u , G 

' Z. pi p + qt q p+ p + qf 

q 


(2b. 4) 


If the system is nearly ferromagnetic, D (0,0) << 1 , one 


should use D rather than D" in such diagrams, i.e. work 
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with skeleton diagrams. It is difficult to calculate Z 

exactly. heverthulcss, since iJ is peaked for small q and to, 

ons can expand the fermion propagator functions in powars 

of ( q ,oj ) . The |q |= 0 , co=0 part or suen a n expans ion gives 

2 2 

the leading temperature depenoonce. i he leading (q or to ) 

corrections give, in the p ar am a g no n liir.it, (r << a ) a 

2 

contribution of relative order t and in the classical 

limit a contribution of relative order and hence are 

small since t<< 1 . This physically amounts to assuming 

that tho coupling between low frequency (to < k n T), long 

\ /3 

wave length (q < /a kp or kf(k^T/£f) ) spin fluctua- 

~y 

tions is of zero range in space and time (constant in q,to 

space). Actually, the time of interaction is 1 and is 

much shorter than typical w The same is true for the 

— 1 - 1 

length scale ( ^ » (2k^) ). One can now compact the 

diagrams 2a-2c into a single term where the value of the 

effectively zero range instantaneous coupling can be 

evaluated through Figs. 2a~2c. It is clear that this 

o 2 

coupling depends on temperature in the scale ( T /Tp ) and 
thus tine temperature dependence is weak. 

Thu contribution of Fig. 2g to the fluctuation 
self energy is 

zdo.O) = - X-b D (q) . 

Pq 


( 2b. 5) 
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In the paramagnetic phase thu- contribution is the same for 
longitudinal as well as transverse spin fluctuations. One 
thus has up to this order, 

JJ { q=D ) = ( 1 - U X( 0) - E X ( 0) ) ~ 1 

= (1- Ux(0) - s J ( o ) ) “ 1 • (2b. 6) 

The zero point part modifies the RPA su sc up t ib il ity at 
1=0, i.o. D T=0 (q=0) = (1 - Ux( o) - ^ =0 (0))"‘. The thermal 

T 

part determines the loading temperature dependence. 

This term is discussed in section zc * and it is shown that 


E-f (0) - T /a Q for - 'o 


t << a 


and 


for a < t << 1 , 

o — 


i-iow consider a higher order diagram, e.g. F ig, 2d, 

Again breaking up the interned spin fluctuation energy 

Integra/ into thermal and zero point parts* one gets 

diagrams schematically shown in Fig* o. The Fig, 3a adds 

to the T =0 self energy. Fig, 3h and 3c have the same T 

dependence as Fig, 2 g 0 iheso are ! vertex ' corrections to 

Fig, 2 a and are of 0(1) 9 due to the inter me d iate couplin g 

(Up ^ 1 ) nature of the system. Clearly they cannot be 
€ F ’ 

calculated reliably. We take the point of view that these 

corrections can all be lumped together into vertices , i.e. 

Fig, 3a~3c can be redrawn as Fig, 3e~3g, Further* since 

these vertices again have a range £p in energy and in 
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( 3e) 



(3h) 


Fig. 3. Higher order diagrams and their temperature 
depe ndences. 
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momentum, the leading order thermal ef tacts are given by 
neglecting their q,a) dependence. Thus the fully dressed 
thermal one spin fluctuation self energy is given in 

Fig, 3:h, where the hatched square represents a constant 
in which all vertex corrections have been included. This 
constant cannot be calculated reliably from first principles. 

The last diagram ; Fig, 3d, involves two thermal spin 

fluctuations. Because of the short range of the fermion 

I 2 

induced coupling this will go roughly as (£y) ? and is 

thus not of the leading order in temperature. However, 
this conclusion is not correct if the two internal fluctua- 
tions are strongly correlated, e.g. Fig„ 2e. Here y for 
external q=G, the internal lines have four momentum and 
— jq ^ . The thermal part of such a diagram can be large if 
the coupling vertex * is sizeable. A similar 

analysis can be done for throe internal spin fluctuation 
leruiSo We shall examine these corrections in sections 
2d and 2b. 

The above discussion illustrates cur approach, which 
can be summarized as follows, 

( i) Because of the intermediate coupling nature of 
the proolem, T =0 properties cannot be calculated accurately; 
therefore, Xj_g and the fluctuation coupling vertices cannot 
be accurately determined. 
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(ii) The temperature dependence of Z arises from 
coupling to thermal spin fluctuations. The temperature 
dependence arising from various spin fluctuation coupling 
processes can be analyzed, according to the number of 
thermal spin fluctuations involved. 

Wc now analyse the temperature dependence of 2 in 
detail, starting with one s.f« diagram and then calculating 
the effects of two and thren internal spin fluctuations. 

2 c . l b pm Fl uctuation 

i 

As discussed earlier this term is given by 

Z 1 ( □ , D) = -X I D( q) . (2c. 1 ) 

P q 

This leads to 

D ( q ) = (1 - U X°( q) + h l DCq*)) -1 5 (2c. 2) 

q.' 

an equation for D(q) which has to be solved self consistently 
for q=0 to extract +he temperature dependence of x(T). In 
order to calculate the frequency and wave vector summation 
in 2 c. 1 5 an explicits depend Li nee of D on q and w is 
rs quireb . We use its functional form in RPA, 

u °ia> = hr— (2c - 3> 

1 - U X ( q, z m ) 

X°(q,i ) is related to the spin susceptibility of a free 
Fermi gas. It is given in terms of the well known Lindhard 
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function F(q,aj). This function has a fairly complicated 
dependence on its arguments. .'Jut for q < < k f and (cu/qv^.) «1 
we have 


X° A (q f w) = (Pe f ) (1 


, 2 iTTYOJ, 

<5 q + } 

qv f 


(2c. 4) 


The si.isl l |q| and to region is obviously the most important. 
In order to calculate X( T ) in various temperature regions 
oi physical interest, we use Lne form (2c. 4) and assume it 
to be valid througii out. Thu error thus introduced is not 
largo. With the approximation (2c. 4) and following the 
discussion in the previous s.ction, we can write the 
equation for D as follows: 

1 


D (q) 


(a ( T ) +6 - 


_iu_Y<i)_ _ 

c i v f” 


( 2c . 5 ) 


where 

a(T) 


= a o + 


,1 


q ' 


(2c. 6) 


and 


a 


( D ( T =0 ) ) 


-1 


This a (and hence a(T)) are related to Stoner enhancement. 


-1 


factors by a” ( 1 ) 


«"!t) 


a(T) is obtained by 


evaluating the integral in (2c. 6) using eqn (2c. 4) with 
the parameters 6 = 1/12 and Y = 1/2 (for free fermi gas) 
This is done below. The energy summation is easily 
carried out. 


1 


I P(q,z m ) = / 


00 h/.i IwD (5»or) 


do) 

If 


r . Bw 


Li 


= ~ / dw 
IT J 


( TTYOJ/qv ) (e^ W -1 ) -1 


. ^ iyw 7 T 

(a +5q ) + ( qVf ) 


(2c. 7) 
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l 

The integral (2c. 7) can be split into two parts,, a zero 
point pare and a thermal part. 

00 1 00 00 e co _1 

/ <J>(to) do) = j <j>(co) dco + 2 / fy[tx>) [e -I) - do). 

—oo Qr _ ] q q 

We have already discussed about the former one. The energy 
integration can be evaluated^ far the thermal part, 
leading to 


8 I E(q, z m ) = C" 1 (in y - (2y) _l - c|»(y) > 

z m 

where 

Y7T 

Q « „.! 

q qv f 8 

y - (a(T) + 6q 2 ) (2 t rk g T C^) 


(2c. 3) 


( 2c . 9 ) 


and ty(y) is the digamma function. Now in order to integrate 
over "q, we need some simpler approximate form of the function, 


4> ( y ) Slny-(2y) _1 ~^(y). 


Let us look at its asymptotic forms. For y»1 


In y 


1__ 

2y 


- 4 >(y) 



•j + _j 

120y 4 2 52y 6 


and for y < 1 


My) = - - + ip(1) + y? ( 2 } - y 2 C(3) + • . * 
In y = y — 1 — 2" ^ y- ^ ) +3- (y- 1 ) - . . . 

I y- 1 I 1 1 ? y ^ 0 


(2c. 10) 


(2c. 10a) 
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w h ic h i mp 1 v 


4>( y) = y - 1 

= 1/2y 


T hus 


<M y ) 


1 _ 

iP 

1 

2y 


i_ .1 

2y + y 

-y < ! - 


y >> 1 
y < 1 


4 >( 1 ) 


(2c . 10b) 


A good rand simple interpolation formula applicable over the 
temperature range of interest is 


In y - 


1_ 

2y 


^(y) = 1 

2y + 12y^ 


(2c. 11 ) 


A comparison between left and right hand side is shown in 
Fig. 4. Using this, the momentum integration can be 
carried out. Before performing explicit numerical calcula- 
tion,, let us first investigate different temperature 
regions of physical interest. From our experience with the 
paramannon theory and the experimental results, we see 
that there are two such regions. 


( i) The Par am aqnon regi me ^ This is defined by the inequa- 
lity y »1 a i.e. 

(a + 6 q 2 ) qv f /2lT 2 y k/r » 1 

which implies T / a Q << 1. Consider the first term in the 

expansion of <f>(y) in eqn. (2c. 10a), and substitute in 
the integral 
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- n l 


l 1 


(2c. 1 2) 


q q M 1 2 y ' 

This integral can he calculated. It is cut off independent 


and the result is 


D(q) 


1 2v, 




/ ! , T \ 

V !n --h * / 


OL 


(2c. 13) 


for x << a Q . This is the T dependence we get in the para- 
magnon Lheory. Note that the calculation is not performed 
self c onsiste ntly „ The leading correction to above result 
can be calculated by taking the next term in the expansion. 


A ( l D ( "q ) ) 
q 


1 

U 


IC - 1 (- i*n) 1 


-*■ q 

q 


120' 4 

y 


(2c. 14) 


— 1 2—4 

The integral turns out to be like J dq q~ (a + Sq ) 

There seerrsto be a logarithmic divergence for small q, but 
this is spurious since we have a lower cut off in this 


limit, i.e. q > q 


(2iry/v f ) (k B T/a Q ) 


Considering this 


we obtain to the lowest order in x 


A(A I D(q)) 

q 


4 


In 


a 


a 


(2c. 15) 


which means that the result 2c . 1 3 is correct if 

2,3, 2,3 

t /a In T /a << 1 . 

o 

Thus the leading correction to the par amagnon result 
is T^/oh In X 2 /a 3 and the approximation converges fast if 
t 2 /a 3 in x 2 /a 3 « 1 . 



40 


( ii) Cl a ssical Curie li ke r egime : Now we investigate the 


other limits y < 1 . I his corresponds to ( 6 q v^,/ 2 2TTykgT) << 1 

cj 1/3 

The fluctuations are cut off at a momentum q = ( — — ) = x 

c k f 

There is no restriction on torpor at u re , except that the 
system is degenerate, i.e.T << 1. Here 


!»(?> -U 

q q 

_ i sl 

IT 2 it 
1 

' ? 

2 7T 


r-1 

u q 2y 

, % 

* k a T > / t 


( 2c . 1 6) 


( a + 6 q ) 


-/£ tan_1 ■’□/!> 


(2c. 17) 


Thus the amplitude of thermal spin fluctuations is propor- 

1 /3 

tionai to (kg T P e „) (Tp £ p' ) . This is the same term 

obtained by many authors^ ’ ^ ^ ^ in their analysis 

of itinerant ferromagnets above T . The term is qualita- 
tively understood as follows. (k^Tp ) is the relative 
J D £p 

volume in energy space of fluctuations with energy less 
than kgT , this cutoff arising from the thermal occupation 


f ac to r ( e^ w - 1 ) 


-1 


The volume in q space of this fluctua- 


tions contributes a further factor ( k D T p ) 


1 /3 


Because of 


8 V 

the size of the coefficient (6(=-.Q8) small) this term 
varies almost linearly with temperature. There is no 
thermal contribution coming from elsewhere, and for a very 
small a(*r) varies as T . This is the desired Curie 
behaviour for the susceptibility. Notice that even though 
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t hs system is highly degenerates ( t << 1 ) , because of the 
spin fluctuations end their interaction a Curie law behaviour 
is obtained and the system uchavss as a collection of 
independent classical spins* Some of the consequences of 
this result will be discussed in Chapter IV. 

2 d * Tw o Spin Fluctuations 

1 8 

It has been suggested by R amakris hrian that the 
contribution of two internal spin fluctuation diagrams to 
a( T ) is negligibly small. This follows from the fact that 
in the classical limit (where only D(q, z =0) is considered 
while performing the energy sum) , for | "q 1 ->0, the contri- 
bution vanishes identically duo to time reversal invariance. 
This is true only in the classical limit but not otherwise. 

Beal --Mo nod et al . have considered this term (they do not 

2 

mention this explicitly) and obtain a contribution x / a , 
same as obtained from one spin fluctuation term. We also 
consider this term. Though in the classical limit this is 
zero but it is nonvanishiny in general. Following 
the discussion in section 2b, wo consider diagrams with 
a momentum correlation between the internal fluctuation 
lines only. These diagrams are represented by Fig. 2e. 

There are four such diagrams, each cont ainingonc transverse 

and one longitudinal fluctuaxion line. Adding all 

these diagrams with proper weight factors and noting that 
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there is no differ cnee in up and down spin in the paramag- 
netic phase, a considerable amount of simplification 
results and wc got, 


E H (0) 


Xvil it o <f>( u) 

2 j?q -3 -a “ 


( 2d . 1 ) 


who rc 


***( q) = 


{ r G k 


:{( 


3 x o+~ 

3B '^=0 


U 2 


(2d. 2) 


1 _ 

3U 


T ho expression in the bracket { )• has boon calculated by 
B e al -Fig nod e t al . 

3 X g+ - 


(■ 


•)„ „ = 2p 


0) 


9B ‘ B =0 " 2 2 

' q v f 


(i+o ((-~) 2 )) + ~ □ (tM 

q v f qv f 


(2d. 3) 

After substituting expressions (2d. 2) and (2d. 3) in (2d.1) 
and carrying out the energy integration, the thermal contri- 


,1 1 


bution to 2 (0) is obtained as 


s n <o)l Th „ 


T hermal 


X, . I / 


dw 


2 ° L J gw , 

q o -1 


^ % {D D ^ Ke ^ 


1 6 P 


2 • v f 2 

L J a.* TT 

e F q o 


2 (a +5q 


q - q 

2. Try to to 2 y 4 


(2d. 4) 


, , rn 2 2 r 2U 2 i? 
^ e F TT^Y V 


2 .2 1 B 


P v f * q 4 


, -2 . 2 ,2 

! TT Y 

2 

} 

+ 2 2 

q v f 


2dt 


, 2nt 
( 0 

1 ) 


; 2 ,22 

{y + t } 


X 


(2d. 5) 
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(a- 6 C; ) qv^ 

Here t = (gw/2iO and y = — -y as before. The 

2 7T'“ kgT 

integration in (2d. 5) is rather complicated, but can be 
transformed to the integral appearing in equation (2c. 7) 
The result is 


.1 


0 I 


k B ’ 


T her mal 


1_ 

sy 


-1 


U) It, , = K /dq — {y^ (In y- ( 2 y ) ' - $ y ) } 


X 2 - a 2 


16U 2 p. 


i .,2 


(2d. 6) 


r 4 ! J ^ 

1 A 


4 2 2 
n y 


■} . 


This integral can be estimated in various temperature 

1 /2 

regimes. In the extreme paramagnon limit x/a « ct it 
2 

gives a y /a contribution, t ho next term is smaller by 
x 2 /a^. In the Curie regime o x^^ term is obtained, as 


expected. 


To summarize we emphasize that the contribution from 
the spin fluctuation correlation (drag) term is in, portent 
and is of the same order as mean fluctuation field term (one 
s.f. term calculated self consistently). 


2e. Three Spin Fluctuation s 

The results obtained in tne previous two sections are 
useful if the fluctuation expansion is convergent, i.c. 
higher order corrections are small. To check this we 
estimate the leading contribution from the third order 
thermal corrections to the self energy. As in the previous 
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section wc consider only those diagrams which have three 
correlated thermal fluctuations. All other diagrams can 
either be absorbed in one or two spin fluctuation vertices 
or they contribute to the order (S^(T))^ or (E^(T))^, 
which is small. The resulting self energy diagram is 
shown in Fig. (2f). The internal fluctuation lines can 
be one transverse and two longitudinal cr all throe trans- 
verse. now there are various ways of joining fluctuation 
lines and the four fermion vertices, e. g. Fig. (5a,b&c). 
Taking into account all the combinatorial factors, we 
can write down the expression Ter the self energy as 


2 in (q) = Xg*-*! D(3i) D ( .9.0 ^ D ( £-J2-| ~S.2 ^ 


(2s. 1 ) 


w it h 




8 


1 511 
8 

1 5U‘ 


-H1JEL2 


_ { ui 1 
2 { 3 t 


;4 } 2 


{ - u }‘ 


Here we have assumed the four fluctuation vertex to be of 
zero range and energy independent. T nis is done by consi- 
dering only q=0 term in the momentum expansion of the 
product of the four Fermion propagators constituting the 
vertex. Writing equation (2e.1) in terms of spectral 
functions. 


„ 1 1 1 / . 
E ( q) 


q l q 2 


P_ ( to i ) 

,q 1 1 

(i)/j z ^ 
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( “2> 


mm z^ +Z 2 


(2s. 2) 



45 



( 5c) 


F ig. 5 . Some 


three correlated spin 


fluctuation diagrams. 
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After performing the frequency summation (over and z^ 

and some algebra, we see that pure thermal contribution to 
III 


Z (0.) is given by 
III 


CO 


E - -(D) = -2j. 3 ^ I_ / H d u . 

R *] ? ^ 2 ^ i 


P 1 (u, ) P 2 (w 2 ) p 3 l “3) 


1 Uj + «2 + W3 


x g^ ) g (to 2 ) g(to 3 ) 


(2e.3) 


iNlow one has three possible temperature contributions from 
the b G se functions g(to), viz. one thermal, two zero point; 
two thermal, one zero point ano all three thermal contri- 
butions. The first and second terms will only modify the 
coefficients of the results obtained in the previous sections. 
The loading characteristic thermal correction comes only 
from the third one. In the classical (thermal) approxima- 
tion with on appropriate cut off (^kgT) i.e. with 

to-| + 0)2 + w 3 ^ kgT , 

r . F ( to ) dto Z 1 k B T F ( (Q) dto 

i c 3w_ 1 5 3to 


and ri ^ 1 ; the energy integral can be estimated to give 

„III._. -- .,2 


z r nk B T 

(0) . -2\,(k T) K ! 


FT P-l ( “-I 5 p 2 ( “ 2 ) p 3 (“3 ) 

nd u . 


q 1 , q 2 


to "I to 2 to 3 


Using 


1 ( iryto/ qv f ) 


P ( q* w ) = n Iw» □ ( q,to ) = ^ 2 , 2 

Ai J,A n 


, . ,TTYto 2 

ia + 6 q ) + (— -) 

qvf 


(2b. 4 ) 

( 2e . 5 ) 
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the energy integration can do performed and the* result is 


Z m < 0) = -2 X-, (k r T) 2 (- T ) 3 


r 

) 


. .2 
4 i ? ^ 2 ^ + ^ 1 


a- 


a + 5 I ^ +"q 2 I 2 
(2s. 6) 


w it h 0 . = 


a 

i 


n k t r* 


0 . 


(a+oq. ) q . j ^ 


T ha tempe-r ature depEindence of S" " ^ comas from a factor 
( kpT ) *“ outside the integral and the angles tan” 0^, Its 
behaviour can now be studied in two different temperature 
re girnes . 


( i) For very low temperatures, T << a 

o 2 


S IH (0) = -2X 3 (k B T) 2 (■") 3 


( 2 r, 




2 

8 h . 


q 2 , 

J c q 1 dq 1 


° (a +<5 q 2 ) 2 0 


q„ 2 1 (n k R I Y TT)' 

J q 2 d q 2 / dt -j-S- 


. 2 % 2 
.(-o q 2 ) 


a + 5 jp + % 


v f q 1 q 2 

(2e.7) 


where we have made the folio', f.ng approximations: 

— 1 —1 
( a) tan - ' = 0^ for two of the angles, since 0^ > tan - 8^ . 

This can at the most lead to an overestimate. 


(b) tan”' 0 _ = 1 for one of the angles. Now (tan 0 ) 


max 


tr /2 = 1 . 57 . Generally, for small T, 63 a',.T << * this is also 
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an overestimate. This simplifies the integration 
considerably. Now the integral over t becomes 


1 


/.; at (a + 6 { qf + q 2 + 2 q q.t)} 


-1 


26 


In 


a+<5( +q 2 ) 

^1^2 a+5(q 1 ~q 2 ) 

which is non singular as a. -> 0 , and a 1 ways positive. The 


singularity at q„ =q is mild and integrable. Therefore 

+ o. 


In 


i-j 


2) 2 ^ 


‘2 
Then 


z In ! 0 i 


A I 2 where I = / 


dq 


ir 


-> (a+Sq 2 ) 2 4<5 1/2 a 3/2 


and m — Constt. T . 


This gives 

E 1 1 1 ( 0 ) = Constt . x 4 /a J , 

2 

which is small by a factor of (t / a) compared to one spin 
fluctuation term in this regime. 


(ii) In the classical Curie regime (os £ T << 1 ) since the 
temperature is higher tan - e factors can be approximated 
as 1. The remaining integral becomes 

E 1 1 1 { 0 ) = -2X, <k E Tf .-S-j-h 

8 IT 


/ d° qi d q 2 


1 


1 


1 


2 2 _i ^ — * 1 2 

a + 6 q-j cs +<5q 0 a + oj qi + q 2 | 


which can be evaluated. 
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The final result is 


(2e. 8) 


£ m ( 0 ) = B T 2 In ( 55 ) , 


(2e.9) 
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in contrast to the linear T. dependence in one spin fluctua- 
tion approximation in this temperature range. 

Thus we see that spin fluctuation correlation correc- 
tion term is smaller, by a factor (x/a)^ in the paramagnan 
regime and by x in the Curi' regime than the mean fluctua- 
tion field term. Therefore, as far as the temperature 
dependence is concerned, the higher order spin fluctuation 
effects are quite small. 

3 

3 . Comparison with Experiments on Liq H e 

There are many metals and alloys which are nearly 

ferromagnetic and our results can be compared with X(T) 

observed in them. Since Liq He' J is a ’clean’ case, i.e. 

there are no band structure complications, we study it in 

2 5 

detail. In their paper II, J. R. Thompson et al. have 

discussed X(T)/C from O.Q3°K to 2.2°K and for ten molar 

3 3 

volumes ranging from 26.5 cm"/mele to 35.5 cm /mole. The 

3 

two extreme cases i.e. highest density (26.5 cm /mole) and 

n 

the lowest (35.5 cnf/raole) ore considered here. These 

correspond to free fermi gas degeneracy temperatures 6.18°K 

and 5.08°K respectively. The Stoner enhancement factor 
5 

a (T) , is a quantity of interest. We transform both the 
observed and calculated numbers to it: 
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a^T)- 1 = X(T)/X p 

= a” 1 ( T ) - 1 

= f (X/L) T° • 

The contribution due to one spin fluctuation is given 
by equations (2c. 5) and (2c. 6) which, for computational 


2 y + 12y 2 


y = , 6 = 1/12, Y= 1/2, T=r/T° and 

it ^ y t r 

q = q/k°. 


purposes, can be written as 

“ 1 (T) = 0(0) + A 1 / ^ dq 

O 

where the 1 b are 1 


A 1 " U_2 


c , v <!> 2 ( up E ) 2 

6 it Y r 


Other functions are defined as 


In writing (3.1) we have made several approximations. 

Firstly, we have assumed that X°(q) can be written as X°(0) 

2 

(1 - Sq ) for values of q of interest to us. For q = 2k ^ 
this leads to a value X°(0J 0.33 whereas the actual value 

is x°(0) V 0.5. For smaller q, the discrepancy is much 
smaller. (We can use full Lindhard function instead of 
the approximate form, but this does not lead to signifi- 
cantly different numbers.) Further, we have approximated 
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the fact that the quartic s.f # coupling goes to zero 
naturally for large q by a very simplu cut off 0 ( q— q ) . 

Now the spin fluctuation vertex form factor can, in principle, 
be calculated in the lowest order, e a g* ? one can calculate 
for a diagram, 



■ • • - — — - — — — 

£-£ 

rapidly for large q. The other spin f luc tu ati0*j coupling 
diagram £. ./X\, 

"O' 






has its vertex going as Y ( 0 E ). This clearly has 

£ £ _P+£ 

a different large q dependence. Further, vertex correc- 
tions are of relative order unity and will alter the q 
dependence. For small q (< k_) and rigorously in the zero 

> V T 

q limit, these two vertex terms coincide. Taking all this 
into consideration, for simplicity, therefore we adopt a 
form which treats all quartic vertex terms identically, 
and simulates the large q decrease by a single sharp cut- 
off. The magnitude of the vertex is treated as an 

adjustable parameter. This is reasonable because the four 
s.f. vertex cannot be calculated accurately. Corrections 

by fitting 


are of order unity. For a given q , we fix 




-1 


-1 


to the low temperature data ( T £ a Q ) where X = X (0) + At 

i . 

ctwM.nL ; ary 


Acc. No. Pk 
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If fluctuations in the same q range determine X at all 
temperatures, and our simple representation of the vertex 
is fairly realistic, we expect a good fit at all tempera- 
tures. os(0) is fixed by X(T=0). i\io n-self co ns iste nt 
results are obtained by substituting a for a(T) in y at all 
temperatures and obtaining a ( T ) from (3.1). Otherwise 
a ( T ) is calculated s elf-co nsis t ently . 


The two s.f. contribution (Eqn. 2d. 6) can be 
written as 


<*(T) = o t (T) 


A 2 t 


/ 


d P 


(1 + 6y) 


(3.2) 


with A 2 = — z°~~2 ^ Up e 

3 it" y F 

Now even with the same simple cutoff model for the two s.f. 


vertex 


V 


, there is one more parameter, i.e. 


the size a 2 the vertex. We fix it as follows. Calcu- 
lations show that the two s.f. term contributes about 
15/1 of the one s.f. term end has the saint temperature 
dependence. So for simplicity, we choose for A 2 "the bare 
value, keeping only A^ as a free parameter. Again self- 
consistent and non-self co nsistent calculations for a(T) 
have been performed. 


The parameters a and A, are given in. Table I. 
r o 1 

The results are shown in Figures 6, T and 8. It is clear 
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that qualitatively there is no difference between one and 
two s.f. contributions if the calculations are performed 
self consistently. Both of these agree with experiment 
over a wide temperature range (x <_0.2). The agreement 
becomes progressively poorer,. The nc e-self c ons istent 
calculation gives good results in the very low temperature 
region only. This is the reason why the results of 
B e al — I'-'i o nod et al do not hold beyond x = 0.2 a^. The non- 
self consistent formula drastically overestimates the 

fluctuation effect for x>0.2 a . In formal terms the 

o 

effect of self consistency is seen as follows. The low 

2 

temperature deviation of ot ( T ) from ot ( □ ) is At /a. That is, 

2 2 
ax AX 

ot(T) _ a(o) = -2 . Clearly for <x( 0 ) » 

a( 0 ) + A? /a 

there is no difference between the self consistent and the 

non-self consistent values. But beyond this narrow range, 

the difference begins to increase. It is also easy to 

see that the leading self consistency correction term - 

A^x 4 /a 3 ( Q ) has t!\_ same a - "' dependence as the three spin 

4 3 

fluctuation term which is also of the form x / a . Thus, 

self consistency effects arc not very important far x<<a Q . 

However, in the classical regime where the corresponding 

2 2 

quantities are ^/a and x In 1 /c i respectively, the 
difference is very large. Thus, it is important to keep 


the calculation self consistent. 




0-0 5 





0.05 
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We notice also that changes only weakly as we go 

from the a = /i to a = ,05 case. This is understandable 
o o 

since A^ depends on Tp or density which does not change 

much. The value of the coupling A^ depends strongly on 

cut off being nearly inversely proportional to it. This 

is not surprising since the fluctuation excitation energy 

goes up only slowly with q (5 small). Thus a considerable 

part of q space is involved, and cutting off a sizeable 

segment of it needs to be compensated by increase in A^ - 

At low temperatures (t << a ) where cut off dec endence of a(t) 

c 

be analytically obtained, the results broadly verify our 
numerical result. 


Now we estimate the higher order corrections due 

III 

to £ . For a Q < x << 1 whore the departure becomes sizeable 

Eqn. (2e.8) and (2e.9) give 


Z 111 = -X 3 U b T) 2 




with 


1 5 U 


2 


3 (2U) 

U p £p 2 
( y) a nd 

24s,- 


6 ( 2k f' 


6 = 


can be written as 


4 

fly In (1/3a®) 0.3) 

0 0 


i/o 9 the correction term 


r 1 1 1 3 ^2 f .4,2,4 2 , s x 

TT (Up e ^_) (j) T In (1/3 a /. 

This gives a positive contribution to a(T) and moreover 

Aa(T) is .02 and .01 at x - . 1 , and .05 and .026 at x = .15 

for high density and low density respectively. This is 

of the same order as A(a_ vn _ - a_, _ _ ). 

EXPT TWO 5.F. 
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I able I 


q 


c 


1 S.F. 


2 S.F. 


q 


c 


1 S.F. 


2 S.F. 




Low density 



v m 

= 3 5.5 chi'/i;io1j 




= 5 . 0 8 0 K 





1 

1 .5 

2 

a o 


. 0870 

. 087 4 

.0876 

A 1 


. 3356 

.2018 

. 1 592 



.0870 

.0874 

.0876 

h 


.3755 

.2266 

.1790 



High 

density 



V 

m 

= 26.5 

cm^/mole 




1! 

— v 

03 

o 

K 




1 

1 .5 

2 

a o 


.041 8 

.0426 

.0430 

A 1 


.3264 

.2120 

.1 592 

a o 


.041 8 

.0426 

.0430 

A, 


.3622 

.2360 

.1956 
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4 . Compa r ison w ith Earlier Work 

There have been several attempts to study x(T) of 

highly p ar am agno ic .netals, alloys and !_;■■, he , but the 

results disagree in many respects. For low temperatures 

13 1 1 

x << a , Deal-l-ionod et al . , K aw ab at a and we find a 
o 

2 

doubly enhanced T variation due to a Contribution from ■ 

spin fluctuations (Eqn. 4.1). These results disagree with 

the singly enhanced T^ term of Stoner's theory. His aw 

3 g 

and Barnea obtain a completely different temperature 

dependence. On the basis of phenomenological Landau theory 

2 

I”! is aw a gets a T In T variation with a singly enhanced 
coefficient, while Barnea in a microscopic Landau theory 
gets a triply enhanced coefficient (Eqn. 4.2). 

X(T) = X(0) [1 - a’ ( - --'- T y -) 2 ] , (4.1) 

o F 

T ^ j 

X(T) = X(0) [1 - a o' In (; r - 1 T -) 3 (4.2) 

fl T 2 a 3 a 0 T F 

r o 

The expression (4.2) exhibits a peak in X(T) at a tempera- 
ture a^Tp/ / g , while expression (4.1) dors not. A peak in 
X(T) is observed in many systems, e.g. Fd, Ot-Hn, L^C^, ^p^°2 ! 

YC 02 , etc. This X(T) maximum is 'explained' as being due 
to a nearly ferromagnetic ferrni system obeying the above 

1 in T 'law'. There are many systems, e.g., IMi^&a, HfZn 2 » 

3 

Liq He , Isli-Rh etc. which can equally well be described as 
nearly ferromagnetic ferrni liquids and which do not show 
any peak. 
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(9b) 


Fig 


9- Ring end ladder diagrams entering in Z1F 
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An expression for AF can be read out from the usual rules and 
the result is 


ff) AF . = 

v ring 2 


y {In (1-u 2 X° + X°~) + l 2 X° + X°-> 

L Q q •-! id 


11 ) £ p 

'V ; r ladd{ 


= l Un (1-U X° + ") + U X° + "> . 

q ■S. 


( 4a. 2 ) 


An expression for c an written down in terms of the 
free energy AF , as 


2 T F 


1 ) 


3 2 AF 


3 a 


3E‘ 


(4a. 3) 


fl=0 


If we differentiate AF twice with respect to B (diagram- 
maticallyl, we get the following set of diagrams (Fig. 10). 
Clearly contributions from one and two, longitudinal and 
transverse internal spin fluctuations, have been considered. 
Ignoring the temperature correction to y^ at low tempera- 

Jr 

tures and considering terms divergent for a 0 , one gets 


and 


((f) i 


3 B 


t AF . ) u _ = + 0 ( T 4 1 n T ) 

2 ring d=0 


(f) (- 


3B 


2 ^^ladder^B=0 


V 

a. 


3 it 


") + 


31T‘ 


TT ‘ 


1 ) a In a ] 
o o 


0 (T In T) 


(4a. 4) 


13. 

+ As shown by Beal-Monod et al . A^F is to be calculated with 
the free energy function evaluated at /a = ^a(T=0). 
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An expression for x ( T ) follows quite clearly. Notice 
that there is no r In T term though there is such a 
term in the specific heat (Eqn. 17, Kef. 13), which is 
associated with the paramagnon contribution to the single 
particle self energy (sec next Chapter). 

2 1 

Recently, Kawabata , using the paramagnon model 

5 1 

diagrams of Ha, !’! c al-l'-io nod and Fredkin, plus the correc- 
tion due to the re normal is ar io n of the chemical potential 
due to spin fluctuations, has also obtained the similar 
results, but with a different numerical coefficient. 

If we also express our diagrams in terms of the 
ferrnion loops, the lowest order diagrams will look identical 
to those considered above. Thus at low temperatures, our 
non-self consistent calculation is identifiable with the 
paramagnon model. It is futile to argue about 

2 2 

the magnitude of coefficient of paramagnon term T /a , 
since it cannot ba determined from first principles. 

4 b . In T term 

Barnea uses the expression (wo ignore etc.) 

= Tr ( < n .> “ < n >) I ( 4b . 1 ) 

B =0 d + " B = Q 

where <n^ > is calculated in presence of a magnetic field 

B . Now 

<n„> = V G 

a p pa 
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and hence 


3 b 

X= l ( 1TB“ ) B = a 

Pa 


(4b *2) 


The Green's function G can be written in terms of the 

Pa 

single particle self energy as 


(v - E. : 

v 1 k O' 


v i - e k -r aE - 2 pa 


( 4b k 3 ) 


T he T In T term comes basically as follows, 


— P — - ( _ n + — -—LEaSjL) 

3B " pa 1 a + 3 B ' 

3G , 9Zjp,3) 

l a gg = l 3 ( - 1 + ° * 

a P pa 


leading to 


x - -2 l g; 


I c s; 


3£ o (p,B) 


(4b, 4) 


Now consider the first term and expand, 


• o2 _ _o3 
i + 2 G Z + . 
P P P 


To .first order in E , this is 


2 l G° 3 £ 

P P P 

2 J( — 1 ) 3 s 

V 1 “ e k P 

o 2 , 

= v ( )} £ 

is 3e k 2 u i - e k P 

a 2 f jf . . 

- I ~^ neZi ^. k) + 



66 


-2 

Since ReZ contains a term In this will obviously 

2 

contribute a T In T term to the susceptibility. And a 
similar term arises from the second part (of Eq« 4b. 4). 

This is tne basis of i3srnea T s claim. 


Vt is now show that th^re is a cancellation between 
the two parts, and hence that there is no T 4 " in T term. 

T hs magnetisation can bs obtained from free energy 
as 


M = < n > 

4 * 


< n > 


3AF 

33 


Consider Fig. 5b for AF. Differentiating this with respect 
to 3 leads to diagrams (Fig. 11), which can be written as 


3AF 
33 “ 


2 2 
.3 E - G Z ) 
PI P4 PV ft 


( 4b. S' 


I {G 

pq 


pt p-i-qt 


ux 


,o+- 


- g 2 .g 


p+ p+q+ <1 _ ux 0+ ~ 


> 


u 2 


3y q + _ 

I nf B — ) - D+ „ 

q 9 1 ~ U X° + 


(4b. 6) 


Differentiating (4b. 5) with respect to 
a term for the susceptibility which is 

3Z 1 

" SB 


I (s 3 z 1 

P l - 


G, 


) 


B=0 


once more, gives 


(4b. 7) 


If we look at the either term in (4b. 7) separately, we 

2 

will see a T' In T terra. However, the differentiation of 
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Fig. 11. Diagrams for’&AF/'fcB. 

! 

! i 
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Eqn. (4b. 6) leads to 


(Ax) 1 - 


o-i — .,2 

y (f x , y 

q 3 B 2 1 - U x 0+ ' 


+ 


i 3 ll + l 

1 313 


) 


2 




U X 


0H~ - 


2 


9 


naither oP which has a T In T term. Since the diagrams 

considered by un reduce to those of Be aJ -Mo nod et al. at 

low temperatures, this conclusion remains valid in our 

theory also. The main point is that a T ^ In T term does not 

appear in a spin conserving approximation for x as is 

obtained by taking the second derivative of free energy in 

a magnetic field (Be al-ihonod et al . ) or an explicitly 

rotationally invariant approximation (as done by us, with 

a vector £ ). The conclusion is true for finite range 

52 

interaction also. Brinkman and Engelsberg have used 

2 2 

U - bq for U(q); the q term will lead to a term higher 
order in temperature but no i In F term. 



CHAPTER III 


5PECIFIC HEAT 


1 . I nt rod uc t io n 

The specific heat is perhaps one of the most important 

thermal properties of condensed matter, since its variation 

with temperature gives first hand information about low 

3 

lying excitations and their spectrum. The T behaviour due to 

Debye phonons, the Dulong-Petit relation and a linear 

relation for the degenerate Fermi gas are a few examples. 

We discuss here the effect of spin and density fluctuations 

3 

on the specific heat of liquid He , again a very 'clean' 

and interesting case. Specific heat in the normal phase of 

liquid Helium-3 exhibits a quite interesting temperature 

dependence. For low temperatures (f << Tp - 5°K) one expects 

a degenerate Fermi liquid like behaviour with temperature 

2 

corrections of order (T/Tp) . The observed behaviour is 

different. The following characteristic features are 

3 1 

observed experimentally: 

( i) at very low temperatures (T < 125 m Q k), C^/RT 
decreases with temperature. The observed behaviour can 
be written in the form C^/RT = a-bT , a function quite alien 
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to Fermi liquid theory. T ne slope b increases with 

2 

density. The curve can also be fitted well with a T In T 

term due to paramagnons, but only over a relatively smaller 

F 3 r 

temperature range" (50 < T < 100 rn k). 


( ii) in the intermediate temperature range C^ 
increases almost linearly upto 0.1 5 0 K and then shows a 
p 1 ate au . 


( i i i ) above 1 0 K 
with a slope of nearly 
value for a free Fermi 


it 



4 ’I h ° 


again start 
R w h ic h is 


gas . 


rising linearly 
half the expected 


All these features are clearly shown in Fig. 12. 


Theoretical calculations can be put mainly into two 

54 

categories. The one formulated by Goldstein is an 
effort to explain the intermediate and ,iigh temperature 
behaviour, Here the total entropy of the system is parti- 
tioned into two components, the spin entropy and that due 
to all other dorr. 33 of freedom. T' — ..y, the definition 
of susceptibility, Goldstein writes an expression for the 
spin entropy as 

S„ ( T ) = -TU-i R 1 n 2 (1.1) 

° X_°(T) 

with X°(T) being the limiting Curie susceptibility and 
X(T) the measured value. The final result for C a (T) is 
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remainder is associated with non spin degrees of freedom. 
Nothing is said about the source of and the temperature 
dependence associated with the non spin degrees of freedom. 
Thus the theory is phenomenological ana merely relates in 
a phenomenological model the thermal and magnetic proper- 
ties of the liquid through an effective u ideal gas magnetic 
degeneration temperature 51 and asserts that the spin inter- 
action of the real liquid is accounted for by the empirical 
ratio X(T ) /X°(T ) „ 

The other category consists of theories trying to 
explain very loi* temperature behaviour, viz., Landau 
theory of Fermi liquid, its finite temperature extension 
and the par amagnen theory. The Landau theory of Fermi 
liquid is based on the description of interacting Fermi 
systems in terms of quasi particles. It is assumed that 
the low lying excitation (quasi-particlc) spectrum of a 
Fermi liquid has a structure similar to that of ao ideal 
Fermi gas. The interaction oetween quasi-particles can 
be described through a self consistent field, acting on 
o quasi particle and produced by the surrounding quasi- 
particles. The energy of the system will no longer be 
equal to the surn of energies of separate quasi particles; 
instead, it is a functional of their distribution function. 
E{ n ( k ) } . For calculation of specific heat one starts 
with the familiar expression for the entropy, viz., 
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I In + (1 - n c ) In (1 


ka 


( 1 . 2 ) 


This formula has a purely combinatorial origin, and its 
applicability to a Fermi liquid is determined by the fact 
that the classification of the quasi-particle levels 
corresponds by hypothesis to the classification of the 
particle levels in an ideal gas. For loo' temperatures 
where very few quasi particles are excited, the compli- 
cated functional eqn. can be simpl if iso and usual proce- 
dure for calculating is followed leading to 


c „ = r k f k fl 2 ' r 


( 1 . 3 ) 


here m* is the effective mass at the Fermi surface; a 
parameter related to the quasi particle interaction func- 
tion. The eqn. (1.3) gives a linear dependence of the sp . 
heat on temp er at Lire , but experimentally C v / T never reaches 
a constant value. 

This peculiar behaviour of the specific heat led to 

3 

the conjecture that terms of the form T In T were present 

in the specific heat. In the microscopic Landau theory, 

using crossing relations for the vertex part and card 

5 6 

identities, Amit and his associates have shown 

3 

that these T In 1 terms are related to non-analytxc 

3 

term of the type |e - y | in the inverso quasi particle 

P 

life time x , which is a general feature of all normal 



74 


Fermi systems and is a consequence of quasi particle inter- 

^ /s 2 

action having a contribution of the form (k.q) 

3 . 

in the limit q 40, This term also lends to (k-k^.) In |k-k^ 
term in the quasi particle energy. The physical process 
responsible for inis term is the repeated scattering of a 
quasi par tide- quasi hole pair. In nearly ferromagnetic 
systems this scattering is particularly J arge and this is 
what produces significantly large results. 

- 4 - 

Using the single particle self energy E(k,£p arising 

10 11 

from emission and absorption of par amagnons, ’ Berk, 
Schrieffer and oth.rs have calculated th>. quasi particle 
entropy from the Fermi distribution n(f$) where e£- includes 

■4" 3 

the self energy £(k,£rO. This also gives a T In T contri- 

K 

5 7 

b ut io n to Riedel has shown, in Pull propagator 

renormalized RP/i, that the above procedure gave only part 

3 

of the T In T term in the specific heat. In addition to 
this part (the 'Fermi' contr ibution ) , there are remaining 
terms which look like the contribution from a system of 
bosons, although they arise from the fermion degrees of 
free dom. 

5 8 

Recently, Pethick and Carneirc have remarked 

3 

that in Riedel's calculation the entropy to order T In T 
is expressed as a sum of dynamical quasi particle contri- 
bution (coming from the poles of the single particle 
propagator) and the Bose contribution coming from the 
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59 

particle hole continuum. However, Balian and de Dominicis 
and Luttinger have established that the entropy is given 
by the quasi particle expression evaluated using statistical 
quasi particle energies. This indicates that Riedel's 
Bose contribution may be regarded equally well as being 
due to the difference between statistical and dynamical 
quasi particle energies. 

These theories explain the low temperature behaviour 
qualitatively but there is controvercy about various origins 
of the T^ In T term obtained. Moreover, there does not 
exist any first principle theory which attempts to explain 
the variation through out the temperature range. In this 
chapter we analyze the effect of spin and density fluctua- 
tions on C^(T). These are the only excitations liquid 
Helium-3 supports at low temperatures. We observe that the 
1 plateau 1 is duo to the 5chottky peak like behaviour of 
the contribution due to the spin fluctuation excitations 

while the slope — — — a t high temperatures is due to 

t f 

suppression of density fluctuation excitations throughout 
the temperature range. It will also be shown that the 
'Bose' and the ’Fermi' terms, both are incorporated in 
the calculations. 

In the next two sections we calculate the contribu- 
tion due to spin and density fluctuations using our 

functional integration procedure and then compare our 

3 

results with the experimental observations on Liq He 
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h Flue tuat ion Interaction model 


In this section we study L^(T) in the fluctuation 

interaction model developed in the previous chapter. We 

-> 

expand Y [f 9 £ ] in even powers of one K odd power 

r 1 

terms are zero because of Furry f s theorem. u The lowest 
order in £ 



gives RRA for the thermuc'y namic potential 

-m-nj 


1. 


= e 


= / (n d £ ) exp {-3 I|£ I 2 ( l-uxt) 

ai "Sr 1 a a S 


ana 


Aft = - * 


3" -> I log ( 1 - UT°) . 
q,z„ ’ q 


m 


This in the low temperature limit reproduces the paramagnon 

3 

results, i.e. a mass enhancement and a T In T term for 
the specific heat. 


The next order involves diagrams wi;h four external 
fluctuation lines, which 

Contributes a term of the form 

q ( q^q 3 3 ? 

+ W 32I3T, 313/5 -( 3 l +32 ^ 



4 * pure f term. 
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The first term gives coupling between four spin fluctua- 
tions while thu second denotes coupling between the density 
and spin fluctuations*. Similarly, the higher order 
diagrams can be urawn and the expression written. Thus 
effectively we have a system of coupled an harmonic oscil- 
lators with a wave vector and frequency dependent coupling. 
Obviously, as such the problem is highly intractable and 
one nas to rely upon some drastic but plausible approxima- 
tions. Let us begin with the fourth order. First of all 
we assume the coupling coefficients K ? s to be and oJ 
independent. This assumption has already been discussed 
earlier and amounts physically to performing an average 
over ! f asL Fermion degrees of freedom ana retaining the 
1 slow 1 5F variables and assuming the SF coupling to be local. 

!m ext consider the coupling with density fluctuations. We average 
over the density 7 degrees of freedom, this gives 


I K b <f*f > C 5 

q q q i q i 

This type of averages <f*f> are discussed in the next 

section. The general result is that because of the density 

fluctuation suppression the average has a temperature 

dependence, far milder than the free Fermi gas term, i.e, 

7 2 

T P ep <<: T * T his will effectively renormalize the 


a b 

The coefficients K and K are written in terms of the 
product of four fermion propagators. 
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’harmonic 1 force constant, but with no extra temperature 
corrections. 


Next consider the quartic term in 5F, Obviously ? the 

standard method is to work in the quasi harmonic approxima- 

4 2 2 

t ion, that is, write x as<x ) x and then estimate the 
error thus introduced,, Doing this the cur ms of the form 


v 2 


k“ <ir„r> ? z * c z , k 2 s <u 


z ! 2 > £ + 5 + 


etc. are obtained. These averages have already been calcu- 
lated in the previous chapter. Thus upic quartic term in 
the quasi harmonic approximation the problem is solved, 
vie have neglected the difference [ <£ > ~ <5 > >]which 
arises from fluctuation correlation effects, e.g., from 
diagrams , where the wavy lines represent 

fluctuation propagators. The distinctivs thermal parts of 
these can be shown to be of nonleading order in tempera- 
ture. We can improve our approximation somewhat by 

putting in the true or observed susceptibility where 
Lt. 


< £ £ > appear. Thus wc work in a quasi- 

q*0 s -q rr 


aver ages 

harmonic approx i i». at io n for the spin fluctuation excitations, 


p 2 

In t he p ar am au net ic p h as e < | £ 1 ■ > = ^ I 5 I ^ and 

we can write the change in thermodynamic potential due 
to spin fluctuations as 

AU= - I l n (1-U X°(q) ^ S' > > 


2.1a) 


-3/e I l n (X p /X(q)) 

q 


9 


(2.1b) 
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After performing the frequency summation one can partition 
the resultant integral into a zero temperature and a 
thermal part. il'Cause of its weak t enm^r ature dependence 
we omit the zero -cemper ature part. its contribution to 
specific heat (linear term) is assumed to be small. 

The thermal part is given by 


Aft 


=~7 l / tan 


-1 


YTTGj/qVf 


1 - U 


X°( q, 0) + 41 D t ( q 1 


■) (■ 


>3co 


-1 


q' 


(2.2) 


• 1 


The denominator in the argument of t an”' function is 

_ -\ 

nothing but the static susceptibility, X(qi)» calculated 

in the one spin fluctuation approximation. (One can 

improve the results, introducing higher order correlation 
- 1 

in X also. While performing the numerical calculations 

we substitute experimental numbers for X(0,0), which 

obviously takes all correlations into account.) Aft 

depends on temperature through the Bose factor and the 

temperature dependence of X t (0,D). The latter have already 

been calculated in the paramagnon and the Curie limit. 

Thus one can get an explicit expression for Aft(T). The 

4 8 

energy integration can be performed exactly, giving 


Aft 


f I [log T(x) - ( x-4) In x 


8 -> 

m q 


x--~ln2Tr] (2.3) 


w he re 


x = 


qv f 


2ir Yk-,J(X (0)/X°) 
B q q 


is* 


qv^ (a ( T ) + 6q ) 


2tt Y k 0 T 


(2.4) 
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Specific heat can now be calculated using the relation 


A C 

' 

g iv i n g 
AC 


-T 


3 k A ft 


31 


k,T 


61 


l 1 


/2 3x 3 x x A , 3 x x 2 a ! / \ \ 

'T Tf + -T> ‘ Mx) + ‘id 

o ! 


2.5) 


where <J> ( x ) = ^(x) + - In x . 


With <J> ( x ) - -j and x given u y equation (2.4) 

2 x + 12x 

one can perform the momentum integration numerically. But 
an analytic estimation of C (T) in various temperature 
regimes is possible. For this one needs a(T) which has 
already been calculated. a(T) enters the expression with 
its various temperature derivatives as 


3 2 

4>( x) — 
3T 


( i ) V er v 1 ow temperatures T <<0l : Here Ct ( T ) is given by 

“ 0 2 

its paramagnen limit form, i.e. a( 1 ) ^ a + h and 

a o 

since x >>1, 4>(x) can be approximated as 


4> ( x ) = -- — “ ' 

1 2x 


with this and observing that a' and a' 1 won’t contribute 
to the leading temperature dependence we find 


k B 


l 


{ 

Ci 


a 


2 

4- U -*■ 

M q 


d> ( x) ( 


3a 

FT 


a ■ 


b 


2 . 6 ) 


AC v 1 

v l 

^ y q (« + 6 q-) 


( 2. 7) 
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which yields AC^^T in (1/a ), the standard paramagnon 
mass enhancement c c nt r ib ut ion „ Next higher order terms 
can be obtained by using a series expansion for 4 >( x ) in 
th is regime, i. a . 

,, V _ _1.„ __J , _1 

x ) — — n + 4 “ g • * ° $ 

12 x 12Gx 252x 
3 

and the first term gives a T In T contribution. The 
coefficient is calculated as below: 


L 


6 krj 

D 


1 

T 


- 


l v, 
q 


1 . j__ x__ 

'3 0 5* u 1 

x H q 


'v - 


1 

1 4 ^ 

( 2tt 

30 

1 ( 2 ir ) 3 

3 

q v 

4tt 

( 2 ^ 2 Tk B ) 3 

0 

T / 
7 

30 

( 2 ir ) 3 v f J 


3 2 . 

q d q 

“7 2TT“ 


; dfl ;*2,3 
q( a + 6 q ) 


( 2 . 7 ) 


2 V'( k r J 


where q = 
^c 


a o w f 


f the low cutoff in this limit. The 


UJ » 

integration is trivial and the most significant term is 

~ — Ur In which leads to a T 3 in T contribution to 

2 a. 


specific heat. This term has been obtained by many authors. 
Later in this section we will show that botn 'cose' and 
'Fermi' terms are incorporated in our calculation. Thus 
we get a T 3 In T term in specific heat but no such term 

in the susceptibility xU). The same result is obtained 

* 

by Beal-honod et al. for paramagnon model. 
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(ii) The intermediate temperature (x <1); Here adj^T/Tp 1 


the ’Curie' behaviour. And <J> ( x ) 'v 


1 x 


s inc e x < 1 . W it h 


these, C (T) can be calculated as below. 

3 ot , i ^ ) 

tures, a(T) aT so that — = Q and = - 

3T d ' 

Using these relations in Lq. (2.6) and <|>(x)wn.~ 


t these tempera- 
a a=D (T) 


T 


2x ’ 


we see 


that we are left wixh a very small term Tor C behaving 


b as ic ai 1 y as ( 1 / i ~ ) ( E q . 2.2) Going b ack to A 0 , , we find 

kgT so th 

specific heat goes to nearly zero. 


that in this regime, Afl ^ k_T so that AC ^ 0. 5o the extra 

D V 


Thus beesuoo of spin fluctuations and their correla- 
tions specific heat has an interesting temperature dependence. 

At low temperature it has a dominant par amagnon like behaviour 

2 

and increases to a maximum and then fails off as 1/T . The 

numerical results are discussed in the next section. 


3 

Now we check whether the T In T term obtained here 
incorporates both the ’Bose' and the ’Fermi’ contribution. 
Let us consider the RFA expression for the thermodynamic 
pote ntial , 


it! = 1 I lln n-UX°fl) + UX°g] 
a 


, l J f~- g(»> lln (1-UX°g) 

q c 


+ UX°ci 1 , 


(3. 9) 
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the contour encircles the real axis from - °° to <» in the 
clockwis- direction , omitting the pole of g ( to ) at o)=0. Here 
g ( to ) is the dost function. The shift of entropy is given by 


AS 


T 


3AQT 

Yt 


l l Un (1-UX°o) ♦ UX°s] 


I / ' 2'ifi g(“) Tx" [In (1-UX°q) + ux°gl (3.10) 

q c 


Carrying out t hu temperature differentiation in the second 
term and performing the contour integration, we obtain for 
this term 

OO 

,2 t du) 


£ u / p(ui) m x + ( “) 

n —oo *• 


3f 3 f 

l ( — Y. -I k - q .: 

f ST " 


3T 


co-e. +e, 

k k - q 


JL tY 


k q 


9(£ k -Y-qx 5 Re X ‘ + (q '^k -^k-q } 


3 P 


3 f 


3T 




3T 


(3.11) 


low wc use thu identity 

3f 

k - Ei - 1 { 


3 f 


g(e u - E 


k-q « 7 L 3T 

(1 - ' (f * 


3 T 


' k-q , l 3T 


9(e k - £ k-q 




(3.12) 
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and rewrite the equation (3.11). Doing contour integral 
in the 1 irst term in (3.10) and combining equations (3.12), 
(3.11) and ( 3 . i ' ) ; 


- l ! V [tan ' 1 

1 + -oc, n 3 T 


UX°( q , to) 


'°l 


-U X l T l q 


03) 


uXr ( q» 0 ) ) (1 - U X°d ( q, 03 ) ) 


( 1 -J X°,j ( q s 0)) ) +(UX°T(q, «)) 


I 


ko 

Yt~ ‘ 


Re yJ e + ) 

a ' K ? ko 


(3.13) 


o 

X I 


and x p are 'the real and imaginary parts of X when 
03 = 0 ) + io . Uorice that there are two thermal contributions 
coming from equation (3.13). One form ? the Boson 

contribution and the second from the thermal part of the 
fermion self energies. Fermion self energies are given by 


T . . , ,2 y / d03 Im X CqfU) 

+ ’ = q * »-= kt + £: k-q t 


< 9 <“> + 1 - f k-q+ ) 


-r 

and (~£ ? 


i 2 1 / ^ 


1 m x + ( q > w ) 


-*■ -co tt 03- e, + e, . 

q k-qt i<+ 


(g(«) + f k _ qt > 


In the paramagnetic phase, we use 


Im X *" ( q » 03 ) = - Im X _i '(q f -w) 


-+ / 


to write 
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^ ^ k > e "j!^ 


= U I 


I'm X ' ( q» w) 


f p* -t- g ( (i) ) 1— f j^~ g^~ 9 ( M ) 

£.->■— c 7»- -*■ + (l) + £ .-*>4 £■£■ •*■ — U) 

b k k + q k k-q 


(3.14) 


This is the self energy contribution coming from the diagram 
(Fig. 11). Thus both the Fermi and Bose terms for AS are con- 


tained in our expression which reduces to (3.9) for low T. 



At low t eriiper atu re s the Fermi liquid supports density 
oscillations also. These oscillations can be generated by 
coupling a scalar probe to the density. The excitation 
spectrum consists of single pair, multipair and collective 
os c ill a t io ns . For low frequencies ( w < qv.^ ) single pair 

excitations share most of the spectral density while for 
the high frequency (w > q v ^) a collective mode occurs as a 
resonance over an almost uniform spectrum of the multi 
pair excitations. F.or a charged Fermi liquid the collec- 
tive mode is plasmon mode while for a neutral Fermi liquid 
(e.g. Liq He 3 ) it is zero sound. All these excitations 
contribute to the free energy and therefore to temperature 
dependence of specific heat. We again work in the func- 
tional integral scheme. To the lowest order, density 
fluctuations can be considered to be independent of 5F . 

The density fluctuation contribution in ’RPA’ is 
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f d c f 

\ q - 

4 TT 


Ilf I (UVX°) 

n ^ J. 


or 


P - k T l In (1 + VX° ) _1 . 

^ q a 


(3.1) 


This is the interaction part, since it goes to zero as 
V-+-0. To this wj should add a 'free* pe.rt. Consider 

|2 

n t — ? it 

3P 0 


= J 


La " l!f q r( p£ P /XD q } 


TT 


i his gives 


<f 2 > 

o 


<P Q P -q > ’ 


X°o 

“pT*' 

^F 


i.e., correct 'free' susceptibility in units of P 

T 

Therefore in this approximation 


fio = r I ln (TT 3 )* 


p e. 


(3.2) 


This gives a T*~ term in which with proper q c adjustment 
reproduces the free Fermi gas value. For V very large (3.1) 
nearly cancels (3.2) and the remaining contribution to the 
specific heat is very small in comparison to free Fermi 
gas value. This can be shown as follow. . The equation 
(3.1) can be rewritten 


n 


l r 

s j 


dw 


gw 


t an ^ { 


V I 


m 


xV<h 


(3.3) 


-1 


1 + VRc 


X°-> («) 


We have just converted the frequency sum to an integral. 

-1 

For low temperatures kgT 't w « e^. and co<qv^.; tan x c an be 
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expanded as power series in its argument and to the lowest 
o.der "the co n ..r iu ut io n to specific ho it is reduced by a 
Factor MV P £ p) , ^ c an be s ee n o n adding (3.2) and (3.3). 

The exact size of this reduction can be obtained as follows. 
Clearly the density susceptibility 

L t . , 1 . 

q-*-C < p cf p ~q > ” 1 + VPep < Pq-P-f > f ree 

and this is simply the (compressibility) ear numerical 

factors. The latter is known, e.g. from sound velocity. 

1 12 

For example, the first sound velocity is ("Irue'^f re e ^ 

^5.0 (at 0.28 at in. ) and M3.0 (at 27 atin.). Thus, one has, 

respectively, (1+VPgp) = 5 and ( 1 +V Pgp ) = 14. Expanding 

the total contribution of density f luccust ic ns to Q as 

-1 

a power series in (VP £ p) , we find that with these value 

TT^T 

of (V P £ p ) , the free fermion contribution is reduced by 

factors of 5 and 13 respectively. Physically the 

3 

reason for thio reduction is clear. I ne He fluid is 
relatively incompressible ; compared to a free fermi gas 
the smallest q (longest wavelen gth* low est energy) density 
fluctuation require excitation energies higher by a factor 
of five to fourteen* Thus these excitations contribute 
only a small fraction of the free fermi gas value to C y , 
and in the first approximation can be ignored throughout 
the degenerate fermi gas range. At intermediate temperatures, 
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( 1 > 1 0 K ) , we arc left with □ nl y t he ' f rc u 1 sp in f luc t u at io n 

f n 

term (th- sf interaction term goes to zero as 1 /T ) 

2 

i.u. 7 r xAi « This is indeed observed,, 

4 . C ornper is on wiv h Ex peri me n t 

be describe here a calculation oT l (T) for liquid 

3 

He y based on chc above theory* Acccrdint to this theory ? 

the specific he et can be written as a sum of three terms; 

- v ( T ) = s - f -(T) + iC v s - f '(T) + C v d * f ”(T). The first is 

o s f 

the free or noninteracting spin fluctuation term C * "(T) 

2 V 

whose value is t C ruL * ( T ) = -jp 'Hr = ( tt‘~ t/ 4 ) . The latter 

T F 

result can be proved by expressing the change of the free 

fermion gas energy with temperature T («T°) in terms of 

density and spin fluctuation excitations with respect to 

the (zero temperature) ground state. doth of these contri- 

s f 

but© equally. Then there is the extra contribution * *(T) 

due to low lying spin fluctuation excitations. This has 

been calculated above (section 3.2, Eq. 2.5) in a quasi- 

harmonic approximation which retains the leading temparature 

dependent terms. The third is the density fluctuation term 

which, ws have argued above (section 3) is smaller than 

2 

Ti T 

the free fermi gas value 1 ^ by a factor of nearly five 
(low density) to thirteen (high density). The separation 
of specific heat into additive contributions from spin and 
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dt i t 
i i ■ ? * • * 
a ’ * ’ * 
vjI * :t‘ 
c , , t , * 1 * 
actii. t 
Lit*' i . 
1 j-n., ,-xr 
z j r j - i «t 


1 actuations is not exact; clearly there is coupl- 

• n tncn.„ This coupling basically renormalizes in 

— urn Independent way fluctuation interaction 

ami propagators* Those cannot be accurately 

{ ~i any jay (see Chapter II, sections 2b s 2c for 

jcCausfi: of the smallness of the density fluctua- 

. , wo ic-lr; it, as a first approximation, to be zero. 

fluctuation interaction term AC^(T) goes to 

i >>aT,V , so that in this range, if our assumption 
r ' 

i .f , 


ruy ard in j is correct* we expect 


h; s - f • i t ) - 


n t 


An t. x ( jn. initio n oi' tho data, c.g. for lot: pressures or 
density ( I ‘ £ 1 ntin.) shows (sec Fig. 17) that this is 


i run. ■ :.i i. rue . 


Ujuv'j T« 1 0 K , the C (T) curve is a straight 


1 in . , ri n d the slope 


C (T) 


= 2.3 - 


4 


Wo thus subtract this 


tort,! from the obo.rved C v (T) a nd exhiuit the difference, 
j_ u l. 13 * f * ( T ) mi units of R as a function of reduced 

* * V 

temperature x= (T/T°). Fig. 13 shows the experimental 
points for low pressure and Fig. 14 shows them for high 
pressure. We note the relatively steep (enhanced) initially 
linear rise. - y (T) peaks at a rather low temperature 
~ n . 02 5) , definitely less by a factor of two or 

pe ak 

three than a Q , the obvious characteristic dimensionless 
spin fluctuation temperature. Further, the peak position 
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v u r v ; , u c h 


r*c: t re] 


u 1 G0 no ^ (within experimental error) shift 
n ohu low and high density limits. Certainly 
linearly to ex . Th v a 1 uc of (AC /H) 

25; the peak is ratmc s.,ia.ll . After the 
esus r at nc r slowly but continuously and 


n c o l» c ii a s i i s 


sbolic in the t ai. 


mion, 


The exp r vision,. Eq. (2.5), is foil no to reproduce the 

breed features or xhe observed C (T) curve correctly in 

v ' J 

e v u r y c tisc, U ri^ input i n f or m a t io n required is the suscep- 
tibiiity inverse, a(T) ? and its first end second deriva- 
tives. i«* take tiiesc from experiment. Obviously , there 
is some inaccuracy in the first derivative, and consi- 
derably more in the second derivative. We have tried to 
do the bust possible by a least squares fit to both the 
observed values and the theoretical numbers. Vie also need 
the ’free 1 fermi temperature Tp. While in actuality the 
spin fluctuation coupling falls off rapidly with increasing 
q y in our zero range model it does not. We thus have to 
suitably cut off the contribution from large q fluctuations. 
This we do by a simple cutoff q , which is the only free 
parameter in the theory. The actual numbers obtained 
depend significantly on the cutoff. This is understandable, 
since the thermodynamic potential is not very strongly 
affected by just the small q fluctuations. For example, 
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n 

c-t In-: t.;’, i\ r.i.tuc,,, In (1 . UX° ). Thus the q-.-O 

R 

f i ui: t out in n'= r nr ..,,Lch 1 -3 X° - a <<l dc contribute signi- 
i ic.-.,nrly ncr.,, t h ± s contribution foils off gently with 

L 1 first, ru.o.'t f ustures which do not depend on the 

cutcM , nrn. trio - . calculation with a ,j articular cutoff 
which giv._r good agreement with experiment. 

f or •; v id w rengo of q c values between 0.5kp to 1 . 5kp 

we f ind that tne initial nearly linear rise, the peaking 

at r: low t unipur atu re , and the slow fall off beyond the 

peak are reproduced. As q^ varies from 0.5kp to 1.5 kp , 

thu initial slope increases, the peak position shifts 

AC 

from =11.015 to 0.03 0, and tho peak value from — 0 . 1 5 

to 0.4(1 (all in the high density case). The fall-off 
beyond the pu uk is slow, and the rate is in all cases 

comparable with experiment. As an example, we show in 

• -r- 3 

Fig. 14, L J *'*(T) with q c = 1.1 kp (hifjh density He , 

T° = 6.1 8 3 K). Wc. used the Lindhard x°^ here ; for this q c , 

there is actually very little difference between this and 

thu small M forni used in calculating X ( T ) . vJe see that 

tho initial incruaso, the peak position, the peak height 

all agree accurately with experiment. The falloff after 

the peak is a little faster than experiment. The overall 

agreement is good. 

s f 

The general experimental features of (T) may 

In the classical regime t>O q , 


be summarized as follows. 
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tht; spin suscl pi 10 il ity is Curie lika, i.e. like that of 

d z^ro temper aturw ferromagnet. As analyzed by f-iurata 
14 

and uonidch, C(T) for itinerant fermion f err o magn at s 

above T rises (due to classical spin fluctuation effects) 
on approaching F^. This is seen, for example, in ZrZn2- 
In our cast, tna rise is smothered by quantum effects 
which reduce the entropy and C , and require (for a 
normal fermi system) that C (T)-*-T as T-»-G. 



CHAPTER IV 


ZERO TEMPERATURE FERRO MAGNET 


In this short chapter, an interlude, we discuss a 

hypothetical (?) case, a ferro magnet with T^nQ-K. This 

can be realized in principle if the pressure on the system 

3 

(Liq He ) is increased enormously (corresponding to Tp 'v 
6.8 3 K), and if .it is prevented from solidifying. 'Here we dis 
cuss only the spin susceptibility. In the second chapter 
we have seen that in the temperature range a o < x «1 , 
the spin susceptibility behaves as 


X, 


X = <x(T) ~ aa Q + T • 

Th'-t is, like a collection of classical spins. We are 
observing a classical behaviour for a n ujnly degenerate 
fernii system. dpin fluctuations have pushed the system 
into a classical regime as far as the magnetic properties 
are concerned. This feature gets revealed very clearly if 
we put = 0 in the expression for susceptibility and 

solve ‘the equation self consistently. As a Q ■* 0, the 
paramagnon regime (x « a Q ) gets narrower and a ’classical' 
behaviour is expected down to zero degree Kelvin. When 


I 
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buch 3 c alculac Lon is performed we get o universal curve 

3 5 ! • calculated susceptibilities for a ^ 0 tend 

a 

c* s y tn r ,i toxic cl 1 y *c o w ar u s t h i s 

This is u very interesting case. Here we are effec- 
tively calculating the susceptibility of s ferromagnet 
with T - 1. T hr quantity of intarest, obviously is the 
critical exponent. The effect of spin fluctuation corre- 
lation on critical behaviour for finite T c fsrromagnets has 

1 8 

.jeon discussed by R amakr ishnan. In the present case, in 
contras t_ to a nun-zero ferromagnet iSFF A gives the 
correct critical behaviour. The reason is the following. 
Suppose for a = 0, a(T) a, x^ . Then the non classical 
( 'uramognon) region t /a << 1 means x << 1 and occurs only 
if X < 1 • This is not possible and so one always has the 
other ( C ur lc -We .ius) region. Here the spin fluctuation 
correlation terra is of the form (Eqn, (2o.9) of Chapter II). 

x in ( ! /a) ^ x In (1/x)^ << x^ • 

If X 'v 1 f the correlation term never becomes more impor- 

tant than the mean fluctuation field term. 

To calculate the index X, notice that theoretically 
for small a 

o 

a = {q c - (f ) 1/2 tan " 1 ^ 

AT T 1 /3 u 4- 

as' a D -> 0 , a(T) « j~ q c . As such q c go-s as T but 
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Fig. 15. X Xn vs. T for a Tr = 0 ferromaanet 
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because of small 5 one can safely put X a T . This is what 
onu uX poets for a purely classical system c 


Vi o have calculated a(T) in a s. IP consistent two 
spin fluctuation appr ox im at io n , The parameter A is chosen 
to bu that of hjujh density case. A linear variation of a 
is obtained (Fig* 15) giving X = 1, for the critical 


exponent. The effective Curie constant is nearly the same 
as that of a free spin 1/2 system of the same density. 

23 

A similar conclusion has been obtained by Hertz 

who has discussed the problem using renormalization group 

methods. The now feature in this problem is the fact that 

1 

because of the Soso factor ( - I) - , the number of thermal 

(classical) fluctuations becomes smaller and smaller as 

T -y 0 (i.u. as the critical point is approached). This 

reduces the phase space for fluctuation correlations. In 

the Ri\l5 analysis, this requires the introduction of a suitably 

scaled 'energy' variable as a degree of freedom additional 

to the three momentum variables q , q and q . In effect, 

x y ^ 

the dimensionality increases by one, and therefore the 
b e h av io ur is me an fioldlike(d=3 + 1=4). We see this , 

Bxplicitly in our procedure of calculating fluctuation I 

correlation corrections per t urb atively . The perturbation ■ 


method converges, the leading correction going to zero 



1 QO 


»nct l.io effect of electronic energy changes 
la r ge. There is no small expansion 


r s . it. r . 


. u me rant ferromagnets << Tp . There are 




... r 33 , 26 66 

i^i-Pt, and i\li-Pd alloys 


u nioignc :iZetion much smaller (in units of jiy' 1 
L ^ 1 number f holes per atom. Gnu can again ask the 

cubic question, ‘.'hat are low lying excitations in this 
system and in o # so they couple to magnetization. In the 
Stoner theory tnese are electron-hole excitations. They 

dec re use the molecular field vjla the thermal occupation 

frb,h7,6R 2 

n;Ct, At low temperatures an unenhanced T term, 

Am -? (T"/Tp^)m :=x^m is obtained. The next step is to 
Cfj no icier the collective longitudinal and transverse spin 
fluctuation excitations. The former are overdamped, 
roonnant excitations and have characteristic energy (for 
low (.,) <vnt‘*ep. Th r » latter are well defined collective 
mooes (pules) with dispersion u) = Dq~ , up to a cutoff 
uj c and then a resonance part. While the importance of 


these excitation for m(T) has been recognized, 


1 5,22 


there exists no systematic study of these effects. In 

2 

this chapter wo show that there is an enhanced T effect 
and a reduced I 37 " 2 term fox Am at low temperatures. 


Near T , a 1 classical’ mean 
c ’ 

\ 


field behaviour is obtained. 



uinne expressing the free energy and its derivatives 


f u nc i> io n al s a 
and 1 * jl 


**"7“ 

£ (Chapters II 


the spin fluctuations 

is an or i i or parameter (Ginzburg Landau) approach, 

its u /• tuns ion to the f erromagnetic case is obvious. We 

separate our the inode E Z which is mac r os c op ic all y 

’oo 

occur ifjc! , .ind . rear integrating over all other s.f. modes 
an: over fermion degrees of freedom, express the free 

energy as a function of (henceforth abbreviated 

as C c ' • For s,r, ali values of £ Q (wea!: ferromagnet) a 
power series expansion of F ( £ ) in which the first few 
terms arc retained is expected to be sufficient. The 

coefficients of £ n in such an expansion contain the 

□ 

effect of the integrated out degrees of freedom, namely 
(coupled) spin fluctuations. The characteristic tempera- 
ture dependences are present therein. However, this 
programme cannot be carried out hers as it stands since 
as we shall see F(5 Q ) depends on ^ in a nonanalytic way. 
In particular, in. .re is a term t j_n c, D arising from 
effect of incoherent longitudinal and transverse spin 
fluctuations, and a ter m * ( T ) 3/2 arising the spin wave 
(magnon) excitations. Because of this, we adopt a 
what different approach where we directly evaluate the 


3F 


equilibrium value of £ Q ? i- e - we eva ' Udtc ^ 

, 8 Cj * __ r (J ) 

equate it to zero to obtain g 0 - s 
can be written in the form, £ 0 - B {< n f> ?Q ~ <n *\ } 


and 

( 1L_) 
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“ '*■'*' d instant and < n a > i s the average number of 

r?i„ ciron^ of s r in a. We include the thermal effect of 

O' in • iuct^atinrio and spin waves to leading order while 

C * 1 c,, Id ting <n > - Further, since F is small, we expand 
0 

in powuro of wherever permissible, and retain only the 

1 ; uii i n y vignn ic rrl terms. The method and results obtained 
tBoc u icc wssed in the subsequent sections. 

next we study the magnetisation m(H,T) in a 

26 L 

gnr.tic field. here, Wohlfarth et ai . have shown, by 

analyzing a GL like model with a 5to ner-Wo hlf art h temper a- 

2 

turn dependence for m(T), that the Arrott plot of m (H,T) 
against ri/m(H,T) should be a straight line. Experimental 
[iJoi.s (u.g, in i*J 1 — P t alloys) show a systematic departure 
for small fl/H and large temperatures. We obtain the 
arrott plot expression in a GL-like scheme retaining s.f. 
effects, and find such a deviation. 

2 . FI uct uat ion 1 1 ite r act ion Theory 

The f err o rna i c phase is characterised by a non- 

zero value for the order parameter (the spontaneous magne- 
tisation). The presence of this long range magnetic order 
enters through the special importance of the time and 
space averaged longitudinal fluctuation field which 

plays the role of a molecular field. We use the fluctua- 
tion interaction theory discussed in Chapter II and modify 
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trt..-;fc Jig effect of (=£ q ) “xactly in the fermion 


t i 


i 'uctuation propagators. The other fields 


r r 


Jl ‘ 1 £•/'* . in will then be treated in the Gaussian 


l1, ' 1 s * iniat jo n . I his effectively amounts to inclusion of 
t i.uc t ant io n e fleet., m fermion self energy to lowest order. 

in treat ine effect of h z exactly we split the 

^oo J r 

/-crimp one ;:t of tad fluctuation field os 


ihuri 


tii on 




C 


- 2 

+ 


q» u 


{■ I / du e- 
2 £ 0 °- u 


( 2.1 ) 


0 + d 


Ii«u |t 


i I 2 = 5 2 

rju oo 


IV® ~ 

« \ I u 


p 

q c 


qu 

5 Z 

oo 


( 2 . 2 ) 


b V 


> \ 


From hern onwar-Js wo denote £-*• by r-1 and 

b qu s qu 

,io w thu expression (tqn. (2a. 11) Chapter II, dropping d.f." 

variables, f 1 for the partition function can be trans- 
cx 

formed to o no involving effect of explicitly. This is 
done by transforming the statistical average with respect 
to an unperturbed density matrix to an average with respect 
to the density matrix inwhich the molecular field (£ q ) 
effects are included in the. Hamiltonian exactly. The 
result is 


Z 

Z 


I 


d J s 

d e. 


-e 


7 . ( £ ) 

"static ^o ; 


/ n — exp {-Y [? a ] 5 /Q } 


a tt 


X 


(2.3) 
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v* he re 

3 

‘•static ( ^q ) = <T U {ex P I" / ° Zs l {u,) du'- ].}> 

o 1 

( 2 . 4 ) 

4 1 

while Y [£ ] c *s calculated in the Heisenberg represen- 
Tj , a _ 

teticn cf (H o -2c‘* S z ). This has the effect of changing 
the single particle energy from ■ e, to E, = (e, + c 2 5 ), 

K K+ K O 

,i.e. there is a spontaneous magnetic field. horeover, 

the fluctuation propagators will also depend on and 

hones the longitudinal and the transverse propagators will 

* 

not be identical. 

We are interested in the free energy difference, 
between the para and ferromagnetic phases. Notice that, 
to the lowest order in temperature the free energy is equal 
to ,(N a-, total no. of particles), plus the ft evaluated 
at y ( T=0) . The proof is easy: let y -r Sy be the chemical 
potential at temp T; then 


F .m i ft (y -+ 5y , T) + ( p +• 6y) N 


unt . " ' 

' y ft( y, T )• + {™ + N) 6 y +, V M 
= ft (y , T) +- : yN, 

■since N PINT is constant in temperature- and p : lays, 

no xclein discussion. Thus one should not consider ths 
;*ein|jBS.ist.urs : .’dspsndsn.i .11. shift whils evaluating F, and 
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*’ w 11 d integration over all variables except*^ is 
« r: rr l . ; ,vn 'irb, a function depending on £ 


i r r ; i 


squetion (2.3) con formally be written 


Si 


i.hiini 


" ; ' e o - G( 5 o ; 

L /? 6 


■ 'S' - 5 r ; * SF e tic U 0 ) + 6F 5 _ f . <?„) 


T n~ sn ini'.’, it fc’t ion condition 


3 u( f ) 

j 


dz 


n 


0 


K 


(2.5a) 


^ + 3 F stt . Uc (C) + BF 5 F U) = ° 


(2.5b) 


yivt-o %(]). Tn obtain 5(T) from (2.5b), one can proceed 

directly, namely develop a systematic approximation for 

p ( r \ ( f .(f) can be calculated exactly) ana then 

IjF'^u static ’ □ 

find the field derivative. We find it easier to write 
down an alternative version of (2.5a) ^ end to examine the 
affect of spin Fluctuations etc. on -cine quantity {<ma>- <n^>} 

occurring in that equation. C 0 occurs i:1 th0 exponent of 
( see e.g., £q. (2.4)! in two places. Firstly, there 

ie the Gaussian term ^ . Secondly, there ie an addi- 
tional term in the single particle energy. - q c Z 

- aba,.). Thus, it is directly clear that 

k4/ ki 


L 



3 t*. ( 4 
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i . 


=> i e 



a k4'~'k4^ 4 


(2. 6a) 


24 ., - Bc >: (<w f> - < n ^ > ) = 0 

o 


(2.6b) 


1 h li r •t ; j. : i : j cte d since 4 q is magnetization order 

t jjr jM.t.ir vyritjj. „ We discuss belovg hoi the equilibrium 
4 ( i ) is c a! c ui r : . ; ri f rom E q . (2.6). 


3 . i > : i y ne t izati i j_n 

i.£! sue P rum Eq. (2.6) that the effect of spin fluc- 
tuations ,.nu Spin waves on 4 q is contained in their effect 
on tin. single particle propagator G p . .-uite generally, 
aim c ,n write 


(3.1) 


J3 + 


r O - 1 
b 

P 


^0 “ 


p± 


In nq. (3.1), I pJ _ is to be calculated in the presence of 
of the molecular" 7ield 5 0 , and of coupling between the 
fern, ion and spin fluctuation vsnablee for exa.npl 

tq. 2a. 6, Chapter II). The usual diagrammatic rules 
(Appendix II) are used. Clearly, can be express 

in terms of spin fluctuation propagators ami fermion 
variables. The simplest approximation for X £ , describing 
the effect of coupling with one spin fluctuation, is 
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(d) 


ri.'j;. 17. The fluctuation propagator D L . 

(a) UFA contribution 1 (b), (c), (cl) one spin 
fluctuation contributions (l> and 3 L) . 


L 



(c) 

Fig. 16 . (a), (b) Fluctuation contribution to single 

particle energies. 

(c) Higher order diagram. 



1 o 


" ’ l£re tns wavy line is a transverse or 

‘ ' i Luc t u at io n p rop eg at a • . We find 


I ^ ' ( P - f 


3 L , “ V -C“J: » G -> 


£' i 


+ e' ^ O 


! rr : f Z. , < n*> 

f 0 


< r > - < n > 

f S, f 


Q £ ' + 

<n^>^ can be written as 


(3.2) 


? /y. 
- -V 0 


f, v.p y } ■ 2 ( b )_ ' ~ v b p G p } (3 - 3a) 


- 2 (^ 3/2 c 3 ill g° 4 : 


i 


+ ft l ff (E - 2 ) 

B ^ ii -E + P~ 

- 2 f~ 5 [r I G° 3 ( Z + Z ) ] . 
/ 3 -u Bp p p+ p- 


(3.3b) 

(3.3c) 


In writing equation (3.3), we have expanded <n^.> ^ - <n^> 

■jb n power soriuu in Z and have kept only the first non- 

vanishing term. This is because the leading temperature 

d e j *u ndenr effect is given by this term. Terms quadratic 

ano higher in Z have temperature defendant ports going as 

higher powers of T (see Chapter II, suction 2b for an 

analysis of this kind of thing), and can be neglected. 

* 

Further* since £ is small we have expanded the regular 


— 1 2 

£ o is macroscopic , of order vH , so that 5 q ~ 0( N) . There 
however need bo no> apprehension that we expana in powers 
of /KL Our free energy is well behaved (proportional to 1\! ) 
and the actual expansion, parameter is a quantity 



1 


1 09 


j t ; rr 


it,: 


a 


~n i r" ) as a power series in £ 

O —■ I + c 2 S Q 

ii “■ 0 


, J ft 


t ici 


i v * j 

£ * 
i 

W I ' iX 


w ;; " i i i < 

" ' * 1 1 r •„.* 

r. i 3 . *1 • > 


f ' ' T 1 , 


; ubic in 


£ o This should be 
s o 

r erroinagnet s . We remark that terms 

J b not be well behaved as functions of 

t» 1 * inn that they are not „ So what we have here 

mi:*.- u i :*jc luuLion amplitude and £ expansion. 

o ' 

can 1.' cn.;v* n that the new temperature dependences 
■'»r’bu f r..r u-ore complicated self energy diagrams 


t ncu u . . i Fig. 16 ) are stronger (higher powers of T ) 

nd hune. neglLgrule for x<<1 (the case considered here). 

Thin is not true close to T as we approach T from below* 

c c 

There t. he spin fluctuation correlation diagrams become 

Lmjior t ant , and rnjr approximations (3.3) and (3.2) for 

< n > - < n, > break down. This critical region is rather 

+ + 

small for weak ferromagnets . A perturbative estimate of its 
width 18 is | x„ - th' 1 ! where x. - (T c /Tp) is the 

dimensionless critical temperature. 


F oat: note ( . . .Lontf . ) 

lO 


of order !\l w , puc! for weak ferromagnets, << 1 . For reference 
we give below the l\i dependence of various quantities occur- 
ring earlier and later in this chapter : — 


T . L, 


£ 'v 
^o 


/U, <n + > - <n,> 'vN, D ’ ( £ ) 


X/ N 


<\, N , 


X; Isi , 


1 

N ! 


U * -j, U . * a 



1 1 0 


/ effect of spin fluctuations 

L,u x 1 * 1 ' zo x i , , at io n ? we substitute the expression 


t.-r 


Eq. (3.3). We m r h e the local spin 


4, “ '" 1,t ll,n c ‘ u ‘ i- «*nij approximation (e.g. writing - - J £ 

* x r : - , ‘ J 'C V S,,' “ ^ I G p 4 ) (I otp, C Q ) J see chapter 
i, .tLChi.cn .-b . nr p disc uss ion ) 7 and find that Eq. (2.6) 


S M - dp e _) c 


3 ll p e F 
o a ^ ~6 ^ 


: p e 

t U (3D L (?q ) + 2D T ( ?q )) = □ 


(3.4) 


In tnrmr of the dimensionless (or fractional) magnetiza- 


1 _ o 

t loti rr< Dx > - < n , > = m = : 

2 1 + /HI 


and the dimensionless quartic 

i n P" Sc 


j | H 0 r 

fluctuation coupling constant u ^ = -(- g—M , Eq. (3.4) can 


bo ‘jritten 


2 

r " u 4 


1) - ™ [3D L (E ) + 2D T (E )]• (3.5) 


In Equations (3.4) and (3.5), the transverse and longitu- 

T L 

dinal fluctuation amplitudes D (? Q ) -nd D U Q ) appear. 
These are 


d l U o ) 


■1 l dL (5 ? 5J 


(3.6) 


and D (5 Q ) 


l 13 (q. 5 n ) 


(3.7) 



11 1 il '•'v'-^on ( 3 , 5 ), the second term on the right 

1 ^-' tuation contribution) is neglected, we 

r ' tk ' ' - >t " nr ' result for the saturation (T=G) magnetize- 

t i ■ ? ' »->. ) orroinegnet , i.e. v 

~ 1 

T '■ • (■••4 ) (3.8) 


!r t! 'jt .ncv- hJ i' art h theory, the temper ature dependence 
r ' ’ i son from the thermal occupation probability of 

■j 

•' uc; Lent-, i.c, "c have, at T/Q, not Up but - U J f £ de 

I 

one thu,> Eq. (j.u) becomes (for T^O) 





( K b T ) 


2 



whnru = 

vJ t ; compare this 
by us later. 



c 


''n° 

TT 


(3.9) 

(3. ID) 


l a nip er ature dependence with that obtained 


In order to compute the temperature dependence of 

m(I) from Eq. (3.5), we need to know D 5 ( ? Q ) • Obviously, 

this cannot be done exactly. I’rom our work earlier, e.g. 

Chapters II and III, we expect that a one spin fluctuation 

T L 

approximation for the self energies E {q 9 £ 0 ) an< ^ £ ^-9* 

should be adequate. We discuss this and more refined 
approximations below (section 4). 



C ^ rom foregoing that cl eterminatxon 

^ ■' - living a complicated self consistency 

1 ‘ r,r ' s on D and iA; ti,:sc in turn depend 
" !! l '* section 4 that, for the weak ferromagnet 

" ' *' ■» problem can be solved partly, and at 

*' ' iCiii i. uv ■ temperature dependences in various 

. <:x lit., ro ru.]i ,r .r,o can be obtained simply. 

1. I 

. *> ;.n.. U 

h nefinitn approxi mat io n scheme for calculation of 
J L ,tnd U i.n the ordered phase is needed. We have already 
c u.ic »;i n r.i.> d then-, for normal phase. There, because of the 
ni.isi uu - . f dpont .intjoua magnetisation and because of the 
i. mt <„t. io r!,i.l invar Lurie a of the Hamiltonian, This 

it; no longer true. The longitudinal excitations are over- 
d un. pod, reaonnnl modes while the transverse are long- 
lived cnJiuctiv! mod© (spin waves) in a small (w/£p, 
q/k|. < in ( 1 ) ) region of phase space ana are resonant 
modus outside this region. We calculate these propagators 
in the fluctuation interaction scheme. 

To begin with, we consider the propagator □*". We 
calculate the St.lf energy in an approximation very 

similar to that used in Chapter II. We retain only the 

one spin fluctuation term, i.e. in E , we keep the 

\ t 

diagrams describing the coupling of the longitudinal s.f. 



1 1 3 


* i 1 - •"•'•I’liL or one transverse spin fluctuation. 

■ r before, to be evaluated selfconsistently . 

" - 1 earlier, the leading thermal effects are 

‘ 1 ■ . “ • u ' i r ’ this approx illation. The correlated two spin 

T - *• r ' i'-i' . pi, per ature dependence is similar and 

■■ '' ' ‘ ci .t; - y a one s.f. term with appropriately 

' 'i' “ ‘ - ft ..onts . In addition, we include the 
CMn:x I: . t Lr, r, r the molecular field £ to leading order. 

-*o 

fir ’v writ - ■ ' ■ . vn the latter term. Tliis is given 
i tiyi til w, ,it ic el ,y by F ig. 17 and has the expression 


mm: 




L’ 1 V / r o r o r o 

2 3 “ (G k+ G _k+g_t + G _kT a _ 


k + q + ; 


U X ( ci) 


3U 4 

U 3 


(4.1 ) 


wtiiirfi 


°k f ( 4- ) 


( V 1 - e k +(--) c 5 q ) 


(4.2) 


The single -article propagators G^, are evaluated exactly 
in pro so r.ce of tho field £ . wow assuming the fluctua- 
tion vertices to be “5 independent, the one s.f. 

* D 

correction term can be written as ( F ig .1 7 b , c , d ) . 


,L 


^ (3D L (£ ) + 2D 1 (£ ) ) 
o o 


(4.3) 


1 , s . f . " U 

W e latqr show how the propag-ator 

au-4-H/ont-. — ^ 


, 1 ^ U , 2 ■ , ^ 4 , o n L o n 1 

D L ( i 3 ) = ( 1 -UX (q) + uf U , (3D +2 ° ^ 


(4.4) 


y * 



J " ' ,,<l 1:1 ‘torr.vB of a relation involving only 

* * ^ ‘ the problem of self consistency 

• ' * i t: ■ . y . 

1 r “ 1 11 -rir,es in approximating u^" similarly. 

1 ‘ : ' 1 r in t ‘ ,~r*: negator has the well known spin wave 

* ^ & 

• * -* * '* ^ J ' Siiiall q and (i) (m our case of q < q c 
u » *'* . ’ * ifi. ,nl need not appear in a given approxima- 

• I . ' r f . r j * I fact, the one spin fluctuation approxima- 

tier. ?';<i - r.ut lead to a pole (i.o. tu (D ) 0 as 

<10 nr/J r -► •] } . show in the next chapter how, using 

U:n .. in entity, one can obtain a spin wave in a scheme 

\* hrrp f L or. t, jrjt inr. effects are included. Here we simply 
uCumu in, it in t,ivu pole region, d"*" ( q , z ) has the standard 

I tl. 


li ( y) 


Per 


„ m Jll „ 

- Do* 


( 4 . 5 ) 


where the spin u av j stiffness D is calculated in the next 

2 

chapter, in the continuum regime (q > q c ^ rnk^., a) > m Ep ) » 
L*{ q) has u broad resonance. Here, ‘w- and q are large enough 
so that m ox m 2 cun be treated perturb atively . (This 
obviously is not to be done in the spin wave or Stoner gap 
regime). Here we write for D (or £ ) an expression 
analogous to that for D L (or £ L ); i.B. we include spin 
fluctuation coupling and coupling to £ 0 both to leading 



S Uir 


s n t 


o f m 


- ii:iut io n smoothly go^s over to the one 

u •' ! tc - r II '’tore £ -+ P Q Since in this 

* 

considered small wilh respect to ii) , q 

T 

- ’ on for D is expectuu to describe the 
• - - - dependent behaviour correctly. We 


v - o - o 
ft L a L 


k^ k+qf 


X ^(q) 


f u A 2 

*<.9> - US 5 0 


(4.6) 


ti'iu.Aiuin is valid for small £ q only), and 


2 


.) . 


(D l ( 5 o ) + 4D T ( ?o )) 


(4.7) 


n h\ c 

Jui) = (I - U X°( q) + yf ? D 2 + u 1 ^ D h^ n ) + 4JJ 1 (5 0 ) ) 


0 ‘ 
(4.8) 


,L 


The equation (3.4) can be used to express D in terms 
nijnct is r.t ion only. 

U L (q) - - U (1° (q) - X°(0)>! 


■1 


2U 


-1 


r '4 r 2 _ 2 ^ i^yoi 1 

1 W 5 o + 5q + qv f * 


(4.9) 


This after integration leads to (Refer Chapter II for 
details) , 



a 


, 2 
tor x << m 5 
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T'"U3 


P**. - , 2 o 

E > -'-'dC Pf 


F m 


T' l/J for t > m 2 


(4. 1 0) 


1 <f -* * positive constants. Since m < ni(T=D), 

* r » means t << m x * This holds at low 

a c 


*\ x •\i ,jr 


* h ’second condition is true only at high x. 


j 

' ’>■< i • nut; t high temperatures, i.C'V (x - t) and 

it ' thiu condition means t >> x - x or x >> x /2. Thus 


u , i 1 - i 


w u 


2 

nx /.ret a T dependence changing over smoothly 


>/.! 


de puntJonco beyond this. 


flu- utTi.-ct of transverse fluctuations has been 
2 2 

.i i-.c udvuJ by Yut.ic.ua also. In the spin wave region he 
.u-L,.in:> ,.s ; _,nd a T 3 ^ 2 term due to spin waves and 

t. hr if in tor ac t in n ’with single particle excitations respec- 
tivuly. Using (4.5) and performing the frequency and q 
su mm, jit io n one u-< tains for the spin wave region, 


d t 


1 


0 


Ep 

1 


,t v s: i 6(» - Dq 2 ) da 

,n(T) ii '"iw T * 

q 


e^ 


iii( T ) (k„T/D) 


3/2 


(4.11) 


(T/T ) 3//2 contribution reduced by a factor m(T). 


-T 


In the resonance region (q > n, Q k f , w> m 0 e f ) D 
is given by Eq. (4.8). It is easily seen that this expression 



* r ^ : that for D l 


Wo therefore find that 


“J t << rn 


* P H; 


/ 3 - 2 

T for T > ,n 


(4.12) 


1 ;,J - ' if - ■ 1 : Mpt r atures the fluctuation contribution 

2 , 2 

'' ; • r i. ’ ... self energy varies as t /£ and hence 

‘■t ’■ i * i c -nc (,3.3c) are not well behaved functions 

2 

■ £ . J me 'i I ’.v., behaviour is valid only for x < m , it 

2 

i- ■ ■ il. lt'M. . T no name term gives a contribution x In £ q 
1 , ■ V'U. energy unt 1 hence the Ginzburg Landau expansion in 
i! iil £ u i nu not hold at very low temperatures. At 
i . i ,’i ; tr mpi r t j I; urn (x - x ), this term does not exist and 


,ansion 


valid as usual. 


:t>, Honolulu this section by summarizing the broad 
f'u _;t. ;ii t*c. In Liu low temperature region the equation (3.4) 


r e ,jd :j 


t 4 - ,r 3/2 

(1 „ a ~y— ) - m (■-”) 

rn m o 

o 


(4.12) 


That is, there iw an enhanced t (just, like in X(T) for 

3 /2 

highly paramagnetic ferrni liquids) and a reduced r 
contribution due to spin fluctuations at low temperatures 
T < T /2 . At high temperatures (near T c ), 


(4.14) 




• "Otica that there is a reduction of 
i ‘Co to spin fluctuations. This point 

' • by 1-ioriyE; rnd Kawabata. 

tl) against T is sketched in Fig. 18 

• , e-.g.j ZrZn 2 ? hi^ba, biPt etc. the 

s i-von, as predicted from (4.11). Further, 
* ' ' flattening out described by (4.14). 



In ; . L . ,, ct ion wo discuss the field dependence of 

'■‘•it,. i : •■■lion. i':dr» is conventionally represented by a 

4, 


■.I wl i V i v ■ i . lor i 
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" i:< i i 


(H,T) against H/l.(li,T), known as 
[i; the realm of molecular field theory 


i . i; ,.i, . .. 1 1 1 1 <■ ' v ■ 1 1 I these are straight lines at I = T Q ; 


l<, 1 


1 .M , 3 1 / hi , 5 

:i ' + 5 ( ~> + 

O O 


(5.1 ) 


Thi ho'jn j',.n to inoicate the onset of ferromagnetism 

:t.'v;ovi..r , recur, tly Wohiforth 2 ^ and others have discussed 
K ( H I' ) for lii-FL alloys. The Arrott olots came out to be 


straight lir...^ over a wide range of temperatures and 
applied fields. There are deviations from linearity for 
high temperatures and low H/M. This behaviour (linearity) 
has been attributed to the correctness of 5toner theory 
for weak itinerant f erromagnets , while the deviation is 
associated with clustering effects. 



u . \ .hlfart'n’s calculation and then 

an effects can be included 0 

; ._, r i or the difference between the free 
'X . Ii.ugnji ic ai LV ordered and disordered 




(H,T) + 7 3(T) i'i 4 1 H , T ) - HM(H.T). (5.2) 


,m ' :s i 


1,.. valid over the whola temperature 
ivromagnets. ninimizing AF leads to the 


X L .1 K lb 


t i-> lor M; 


O.U) + ■■{>.) il(H,T)) M(H f T) = M 


( 5.3) 


i i .. i 

Ad) 


• ■ 12 X (n,T)l 


-1 


U * '• ■ 


■ !: ( n - (;:*h 2 (o,t) x(o,t))‘ 


jjhu x -x 


X(ri,T) - [3..(H,T)/8H] H(T « 

using the 5toner result 

2 „ 2 


■‘•he high field suscep- 


t ib i 1 i t 

M 2 (G,T) = i-riu » 

end X(O.T) = X(0,0) 


2 / _ i- 2 (0,0) 


-2 ,-r 2,-1 


one gets 


2 X^H 


( 5.4) 
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~ hiS is tf1e Arrott Plot expression extended to a wider 
tsr.-eratjre range. The slope of straight lines ■H 2 (H,T) 
o *j 3.1 n;» t ft/:-. :.a independent of temperature. 

Tw consider the effect of spin fluctuations on 
( 5,4) t we incxuds the magnetic field effects in the 

is x /.re s&iafi (2,3) for the partition function. This will 

yPo 

od;:, a torn *—'»<* H(=h) to the molecular field"? . Therefore, 
z ~ o 

the analysis of the section 2 can be applied directly, 
leading to 


2 % -& C * ‘<"<>1 +h - < n *- V +h ) = 

* n 


c- £ 


This is the equation f.or"^ o (h,T) (or m(h,T)). Including 
the effect of fluctuation self energy to the lowest order 
and carrying out the calculation similar to section 3 
and 4, we find 

U 

. (T) + 

T 


m 


- 'Uf^ (T) + (20 . + 3D ) + 


+ 2U 4 (m+h) 3 J =0 


, (5.6) 


Retaining, only terrn linear in, In and using the expression 
(3,5) for m(0,T) we get 

h(1 + (! 4 m 2 ( 0 ,T)) „ 

■ «»( h,T)^ r m{ 0,T) > r . ^ j iT H T TT^ , 

T his expression will obviously give a different. m(.h,T) . 

.Now .there:’, * n .*.*? d •“ P ?. C 7 - 

" * « , J/l'L ; 'itev't '".iambi's kiuie s't'hie^gryesr an - e *f**^io^ V 
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.4/3 

varies as T 9 the parameters fixed ror low t em per ot ure 

points will not give agreement in the high temperature 

2 

regime ? and hence the deviation. Since rn decreases 
with temperature, the slope (m(H,T)^ - rn( 0 f T ) ^ ) /( H/M( H , T ) ) 
will decrease with increasing temperature, in agreement 
with the experiment (Fig. 19). 



CHAPTER VI 


SPIN WAVE STIFFNESS 

1 . I ntrpduc t ion 

Spin waves exist as well defined, gapless, collective 
modes in isotropic f err omagnets . Their existence can be 
established in exther of the models of ferromagnetism, i.e. 
Heisenberg model or an itinerant electron model. In fact, 
they are associated with the breaking of rotational symmetry 

m 

at the second order phase transition point, T . The energy 
of this mode for small wave vector q is given by 

_ 2 

0 ) = Dq (1.1) 

The coefficient D is known as spin wave stiffness. It is 
related to exchange integrals etc. As the temperature of 

the system is increased above 0°K, the energy of the mode 
keeps on decreasing. Above the transition point there are 
no long-wavelength spin waves. Since we have seen earlier 
that the temperature dependence of physical properties in 
weak ferromagnets is dominated by spin fluctuation effects, 
it is obviously relevant here to discuss JD ( T ) and effects 
of spin fluctuations on it, in an itinerant electron model 
for weak ferromagnets. 
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We first discuss earlier theories, for orientation, 


In the RP A < 


the transverse spin susceptibility 


X (q*oa) xs calculated in the ladoer a,-, rc-imation. It 

develops a pole when the zero range repulsive interaction 

l) attains the value implied by the Stoner criterion U p_ =1 . 

e F 

For larger values of U there is an exchange splitting. 

With vhe electron propagators appropriate to this polarised 
st ate 

-1 


k + 


(v 


1 ~ e k — 4/2 > 


( 1 . 2 ) 


the susceptibility becomes, 


X 




X ( Q.to ,A ) 


UX ("q ,co ,h ) 


(1.3) 


Here X is the generalization of the F auli susceptibility; 


X j I b j< + qt * 

The character of the singularity in Xiq> u ) for smell |q| 
and u is the i determined by the Golasione theorem. 
Specifically, when the ground state spontaneously breaks 
the symmetry of the Hamiltonian via a non.-zera A , that 
symmetry must be realized by a Goldstone Boson, whose 
energy goes to zero in the limit of infinite wavelength. 
That is, a uniform rotation of all spins costs no energy. 
This implies that X + ""(q,a>) has a pole at zero q and w. 


We have 
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1 = UX°(O j Q 5 A) = - - 

a e 


( G 


k f 


G kT } 


u 

I rn = 


where m is the magnetization. This equation determines A 
which is just the Hartree-Fock exchange splitting. The 
Goldstone coson is s spin wave and we obtain its energy to 
as a function of q by solving for the poles in x(q>w’) 


1 - UX°( q.to. A) = 0 


•(1.4) 


For long wavelengths the result is 


-(5) -i III ( 


+ n ,.x) (q- 


3 ,2 


2 v "kt T ,, k + / 3l< 


( n |<;t " r 'l<4' ^ 


U*>n 


( q 


3e 


k > 2 


3 k 


r i 


(1.5) 


However, if the q value of the spin wave is large enough, 
greater than certain q m 9 this pole will merge with the 
continuum, i*e. spin wave decays into an electron-hole 
pair and is a resonance in the magnetic response. The 
details of spectrum of these transverse excitations are 
shown in ’rig, 2{j . 


For a cubic lattice Eqn„ (1*5) gives 


D ui { T<*lAt v 2 e _ 

3 k 2A k k 


kt ' - k ~- (V k E k ) 2 > (1.6) 


Inserting in this aquation the values of one electron 
energies, one can investigate the dependence of D on band 
structure, electron number, exchange splitting, magnetiza- 
tion etc. For weak ferromagnets (m << 1 ) , by expanding D 




A/k F 


Fig. 20- Electron - hole excitation spectrum (transverse) far 

an itinerant ferromagnet. p 
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in powers of magnetization, a linear relation between 0(1)" 
and m(T) has beer, obt ained . J This relation is analogous 
to the dependence deduced from the Heisenberg model for 
insulators. In general, D goes to zero simultaneously 
with the exchange splitting A at T . Using the same eqn. 

17 

I'iathon and ;.ohlferth have shown that 

D = D q + D^T 2 , (1.7) 

and have calculated for tight binding bands. For this they 
have used the Stcner expression for m(T). 

9 

Recently Young has improved the HR A expression for 

D taking into account the correlation effects in the T- 

2 

matrix approximation. A T term for D(T) has been isolated 
by using the Sommerfeld expansion of various temperature 
dependent functions (e.g. the single particle occupation 
number, the average t-matrix etc.). He has evaluated 
and for a f 0 c 0 c 0 lattice with first and second neighbour 

i 

inter ac t io ns . 

2 

Thus in R'r,\ the stiffness has basically a T depen- 
dence because of the thermal single particle excitations, 

2 

but the coefficient (^p e ,_) is too small* by orders of 

r 

magnitude, compared to what is observed in, say* Mi. 

2 7 

Izuyama and Kubo have discussed the temperature 
dependence of D on general phenomenological considera- 
tions, as follows. 



To lowest order the Hamiltonian of interacting spin 


waves can be written as 


i? k 


«, - I 


kk 


n k n k' C kk> 


( 1 - 8 ) 


is the number of spin waves with energy 0.^ and 
is the quadratic coupling. In the mean field approximation 
one decouples t.ie second term, getting 


H = l (£2 Tc - 2 l C ffi' <n R" > ^ S \ E * ° 


(1.9) 


In the interaction term in the only dependence on k 

arises from the coefficient C, , , . C must qo to zero as 

!<k ’ v 

k -*■ 0 , otherwise E, would possess an energy gap at k=0. 

2 

Because ws expect k at small k and since C^, , rs 

2 -> 2 

symmetric in k and k' , C^, a. a k k' . On taking 12 = Dk 
one has 


E k = k- {D - 2a l k' 2 < n^ ,>} 
= k 2 { D - a' ( k & T /D) 5/2 > 


( 1 . 10 ) 


to first order in a . Thus spin wava interaction gives 

5 /2 1 

rise to a T " term in D(T). This dependence is the same [ 

73 : 

as in Heisenberg model. Obviously so, since the assumptions 

made are exactly/ valid for either model. 


The effect of interaction between spin waves and the 
thermal electrons has been considered by writing the Hamil- 


tonian as 



n n C 
ft k ' Rk> 


H el,s.w ” n |< % ~ J ( 

+ l " h' Er ^' 5n k’ 6n R' 

~^^Rk ?n f<^ n k’* i 1 1 ) 

k 

5 is the number operator for the thermally excited elec- 
trons and C is the electron-spin wa«e coupling. Again in 
the m. f . a. , 

E i? =■ %>- 2 c bb” ‘V* - c nb' <{n k’ > ■ 

2 

The last term will obviously give a T contribution to D( T). 
He nc e 

D(T) = D q - D^ 2 _ D 2 T 5/2 . (1.12) 

The above qualitative prediction is supported by 

7475 

neutron diffraction results on F'e and ;’m i . ’ But there 

is an ambiguity about the exact power law obeyed by D(T). 

29 

For example ^ the spin wave resonance results suggest; 
say, for Fe above 100 D K, a T y depenaence for 1 - Oj/Li . 


5 imil ar ly 

for di--Fe 

films 

the 

s wr 

results agree with 

D = 

D ( 1 - a 

o 

t 5/2 } 

for 

T < 

80 D K 

and 1) = 

D 1 (1 - a 

I 3 /2) 

for 

T > 

90° K . 


Moreover, the observed change with temperature is much 
larger than that predicted by theory. The conclusion is 
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that there is a need for investigation of D(T) beyond mean 
field (KPa) theories. 

The only inter ac t io ns t ake n int^ Recount in RPA are 
the scattering of single electrons against the Fermi sea 
of electrons of opposite spin (the molecular field of this 
problem) and the resonant scattering of electron and hole 
of opposite spins responsible for the spin wave pole in X. 
Interactions between the single particle excitations and 
the spin fluctuations are ignored. As elsewhere, these 
should have important effects here also. The problem is 
to incorporate these effects in a spin conserving manner. 

The standard way of doing this is via a bard identity. 

Recently there has been progress in this direction due to 

76 

Hertz and Edwards. Using W ard identify they have discussed 
the effect of elec tron-magno n interaction on the structure 
of single particle self energy, spectral functions and 
spin wave stiffness constant at zero temperatures. be 
modify and c ppiy their analysis to the problem at hand. 

2 • F orma l Re sults a la 'riert z-Edw a rds 

The spin wave stiffness has a simple relation^ with 
Y(k+, J<-q_f) , ths irreducible part of the spin flip current 
vertex T ( _k 4- , _!<-£+) for coupling to an external spin flip 
field : 

m DqJ * C(q) - ’aJ l ■ k -+> G k + G kt (7 lf E t>S \ \ 

( 2 . 1 ) | 

1 

| 

• i. 
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where 


C(q) 


'k+q 


= q .q . y 
2 i'j ii L 


e h <n l<+ " n '“' 

3 2 e, 


k+ qf 


k 


a k i 3k <n kf ‘ n k4- > 


( 2 . 2 ) 


to order |q"| . Thus the knowledge of y and the Green's 

4 - 

function suffices to determine D e We have to calculate Y 
and G’s in an approximation scheme. The scheme should be 
spin conserving so that the spin wave have correct 
dispersion at low momentum. Edwards and Hertz^ show expli- 
citly that the satisfaction of Ward identity implies that 

2 2 

spin waves will have a) =Dq and not oj = A + Dq . The Ward 
identity is a relation between the vertex T and the propaga- 
tor G . An arbitrary approximation for G (i.e. the single 
particle self energy 2) and an equally arbitrary approxi- 
mation for T will not satisfy the Ward identity. There 
are standard prescriptions For choosing the vertex T given 
a G or Z such that WI is satisfied. The choice of 2 is 
to be motivated physically. 


is 


The Ward identity as obtained by hertz and Edwards 

- 4 * - 4 - 

q y r y ( j< 4- , k.-.q+) s q o r o (k+, ]<-.£+) ” q * T + » k-S+) 

= G" 1 (k) - G" 1 (k-q) . (2.3) 


We analyze the equation (2.3) closely and try to get a 
simple criterion which follow as a consequence. We first 



define the various terms involved. T is the vertex for 

o 

coupling to an external spin flip field. It is related to 
the transverse susceptibility by 

X(jl) = ^ | r o ^ — —"'•2+ ) 6^(]e) ^ J< — id ) (2.4) 

F o is reducible. Consider Y ( J< « ? the irreducible part 

of F - Figure 21 shows the relation oetween T and Y * 
o o 

If we divide any reducible diagram for T at its right most 
part of reduc ib ilit y 9 everything to the left of that point 
is a contribution to full X(q) and everything to the right 
a contribution to Y(J< ? J<~q)« 


r 0 (k.+» j<-q+) = X (k.4- , „k-£t) + u x(.q) Y(k.t> J<-^+) (2.5) 


A similar analysis can be applied to the diagrammatic struc- 


ture of the perturbation series for f , giving 


j<-£t ) = y (Jk + » _k-_q+ ) + U Xjg(ci.) Y(k+, k-qt) (2.6) 


where y is the irreducible part of f and Xj^ S P 

density current susceptibility: 


X ,c = <T ( Jl ( u) 5 Z ( u f ) ) > 
Jb q jq 


which satisfies 


Wl = B l r(J<+, J<-a+) S _ K+ G K q+ 


Xj5 an "^ X » which from 


We need one more relation, between 
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Fif 


21 


The structure of the reducible vertex T~. 


I 1 ^ 

( a ) 


f 


t 


f 'T 

(b) 




* 



(<*) 


Fig. 22. The lowest order class of self energy corrections 
due to (a) transverse fluctuations (Ladder series) 

(b) longitudinal fluctuations (Ring series) 

(c) Hartree term. 
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the equation of i, lotion of X(q,u-u i ) s is 

q . Xj 5 ( q) = q Q X(q) + m (2.7) 

Combining equations (2.5), (2.6) and (. 2 . 7 ) with (2.3), we 

find 

G k4- " G kl q 4, = ^ q 0 " Y( j<t , k.-S + ) •• q • Y (JS^ > 

which at q=0, q ( =co) = 0 reads 

o 

Y(k+, k+) =1 - (L i+ - E k _ + ) /A (2.6) 

Here we have used 



The equation (2.8) gives a relation between 2 and 
the scalar y obeying the W ard identity. Once we get 
diagrams for y corresponding to an approximation for E 
we can construct diagrams for the vector vertex y and 
hence calculate u(T). We use the equation (2.8) to check 
explicitly whether our approximation scheme satisfies the 
Ward identity. 

3. The 5 elf Energy a nd Scalar Vertex 

Vile have first to choose an electron self energy 
function Z . Since the vertex is obtained by differentia- 
tion of particle or hole lines in diagrams for E one 
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necessarily requires an explicit diaqrainmat ic form for E . 
This somewhat limits the (temperature) range of our results ? 
since the self consistency present in our earlier results 
(for l'\ etc.) will no longer be true. A/e shall 

ootain results fur D in an approximation where the leading 
fluctuation effects at low temperatures are correctly 
incl u ded „ 

So far we have been working with the Hamiltonian 
involving electrons plus interacting fluctuation fields. 
Since we have to work with electronic degrees of freedom 
explicitly in order to formulate £ , y etc., the logical 

thing to do is to integrate out the fluctuation field to 
get back the interacting electron system. This gives 

H = T + V , where V = - jU S.5 

Except for a one body term, the potential energy can be 
reexpressed in terms of the well known Hubbard model, 


H = T + U 


n . . n . . 
x+ X* 


(3.1 ) 


The spin fluctuation effects are contained in the 

propagators D T and U L (or x + “ and x* 2 ) . For a weak 

f err omagne t , the farmer has a low lying resonance (energy 

2 

m £p) while the latter has a spin wave pole, again low- 
lying because of small stiffness. We therefore use for £ 
an approximation where the electron couples to X ZZ and X + ~. 



w aves w it h to 


find 


= Dq as already discussed. 

Next considering the diagrams Fig. 2ra and b, we 


E = — TU 2 (X + ”G, + X" X )+U<n4-> (3.3 a) 

k\ B £ X q i<-qi £ J<~.q + 


and 


, 


k j, = -- I U 2 (X + G. + X + n G k _ ,)+ U <nf > (3.3b) 

j3 “ q j< - q f £ J<~£! y 


where 


\a 


X 5 = i - u x 0+ - 


x r 


04 * 


4 * 


x s ' , :~ii‘ f x° + x°- 

3. 3. 


x 0 + -(q) - - J \ V Vat 

and X 0± ( 3> = " 8 5, & ko S k + qo ’ 

The vertex Y is now obtained by inserting h (wherever 
permissible) once in an internal bare line in the diagrams 
for the self energies we have just calculated* ^e get a 
set of diagrams represented in Fig. b«f„ where the 

hatched boxes represent one of the terms of the RPA ring (L) 
or the ladder (T) series and in f the double wavy line 
represents the same ring series with the first term being 

ns v /y ^ 





These diagrams represent the following expression 


for y(Jk4-, k.-q+) 


? I £, , G. , G (• 

8 2 Jk f q T “ + «• V il + £ + 


X ~. • q-q 1 > r 
. i , y 0+ - ^ V+S' + 

V 


^ ^ , G i< ' T G J< ' - q t G j< +£ ' t G J < 1 - q- q It- 


u 


(i -ux°r-)d 


h l 


G G 


■u 3 x+, 


6 k'ja' -k+q'f k’+q^ {1 _ ux ° + - j} 


-a-a.’ 


e; I G G G. , G. 

3 k t J< 4- j< -£f M l Jk-q- q_- + 


U 


(1 - u x° + ~ ,)(1 - u 2 x°T x°7) 

-q-q ' q ' ^q ' 


j_ v 

6 q, k-q ,+ + 


-u 3 x°t x^T 


(1 - u 2 x°t x°7) 


1 l 


We now check whether above set of diagrams for E 
and 1 Y satisfy Eqn. (2.8). This can be done after some 
simplification and some algebra. The algebra is straight- 
forward, the main simplification is in writing the product 
of three Green’s functions as the difference between the 
free particle like susceptibilities. For example, consider 
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3 k» G -'"a f + G is' + 


which for c[ = □ , reduces to 

i i i - 1 + i-' f 


« i y 

3 


k’ Um + V+‘ Jk'+ ‘A' +J 3' + 

< 4 * & k , + - G + ) 


K ! 


1 

IT 


(x'; 


m 


m 


G + — v 

X,,t ). 

lL 


J< * +q ? + 


This is correct to the lowest order i » i j, , which is what 

\< CT ? 

we consider. After such simplifications, adding all the 
diagrams for Y and E separately we easily get the desired 
result 


Y( k + , k t) = 1 


Z, + - E, , 

Jsf kt 


This confirms the correctness of our diagrammatic 
expression for Y starting from the approximation we have 

made for £ . 

a 


4. T he S ector Ver tex and D(T) 

We consider a set of diagrams for Y exactly similar 
to that for y . These diagrams can be classified physi- 
cally into three categories; namely, the self energy 

correction to the RPA propagators, magnon drag and the 
vertex correction terms. In RPA, Y was unity and the 





propagators used were in HF A . In the first category the 
lowest order spin fluctuation self energy corrections are 
included (Fig. 24b, c). The magno n drag term is represented 
by diagrams b-e in Fig, 23. Physically, these diagrams 
represent inclusion of effect cf fluctuation interaction 
or the ' magno n drag’. These diagrams are similar to two 
spin fluctuation diagrams considered earlier. These are 
also included in the phenomenological theory of Izuyama 
and Kubo . A typical diagram may be represented by 
Fig. 24d. The next set represents vertex or electron-hole 
wave function modification by longitudinal spin fluctua- 
tions ( F ig. 24e) . 

Now we consider the Hertz-Edwards equation (Eqn. 
2.1) for stiffness. In RPA, for weak fe rromagnets , the 
stiffness is proportional to magnetization. Here also, 
we expect similar results to hold except the fluctuation 
corrections. That is 

m Dq^ = A(T) q^ (4.1) 

to lowest order in m. This can be proved as follows; 
Consider relations 


mD q = C ( q ) 


[ lim lim 

^ * w -K3 J q j-*- 0 J 


( q s 03) 3 .q 


(4.2) 


and 


^ y y . ^ hp «*)»- |„y 

w[oix(q»w) + m] = - C ( q) + q . Xjj(q» w ) • R 


(4.3) 



from HE. Since 


.Tim lim 
W -+-0 | q }0 


2 


9 


Ml 


m*+-0 


kial 

= 2 

r~\ 


lim lim 
to-* 0 j qj-s-0 


1 im 
nr^O 


X 

JJ 


( q,to ) 


C ( qj. 
2 
q 


c(_q) 

2 " 

q 


2 2 
03 — g , 

03 


1 im jq"|-»-0 
lim 0 ) -* 0 


__ n 

= u 

2 

Thus mDq does nox have any term independent of m « The 

2 

lowest order term should be of order m # Therefore we expand 

each term on R!!S of Eqn. (2.1) and examine the thermal 

2 

contribution to the coefficient of rn term 0 


We now analyze D(T) considering various contributions 


to y 


a * 


Bare RPA 


This corresponds to Y = 1 and hence 


y 


a 


a * 


The propagators G, are the Hartree F oc k one, 

K (J* 


ka 


, A.-1 

= (v, — e_> + a — ) 


T his gives 


q . f - I y Q« + » 

k ~~ 

= q 0 q B "b" I k o k e { 

Expanding in powers of A , 

2 

A , we find 






_1_ 1 

v - + A?2 v x - ej* - A/2 

the lowest term is of order 


} 
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m( T ) D 1 (.IP A) = 


36 


(4,4) 


This is the well known result for low magnetization. 


Sel f e rmer cp/ t corr ec tions Here we dress the Hf propa- 
gators with the transverse and longitudinal fluctuation 
lines. This is represented in F ig . 24 c & c and contributes 


q a q (3 3 ? k « k B G k + G ki tZ kf b k + + E kT b k^ ] 


to lowest order in £ . We have already discussed the 

approximation scheme for E’s. 5ubstiluting those expres- 

2 

sions, we find to order A 


q a q B a2{ 3 L l k a k $ E ° 6} { l f""(q) + X + (q) + X“(q)> 


+ qci 


2 

L k ct k 3 b kh {“- l ( X + ( q) - X~ ( q ) } 


B % 


C 1 


(4. 5) 


3 /2 

The integrals will give a m(T) (T/D) ' contribution from 
X ( q) (due to spin wave pole) and p t /m term from the 
longitud inal fluctuations. This has already been discussed 
in detail in Chapter V . 

c. Vertex cor rec tion ; This is given in Fig. 24 f and 
contributes 


k , l. *„♦ * k + <w + 


B k 


k+jg 1 t _k + q 1 + 


x { 




- U 2 x°t x°7 


■} £j? • q 


i 
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We assume _q' =0 in comparison to the electron momentum 
energy _k wherever possible so that the integration can be 
divided into two parts one involving sum over the product 
of four fermion propagators and the other involving sum 
over fluctuation term. The first p t then again can be 
expanded in powers of giving 


( mDq 2 ) 


v srtc x 


1 a 2 2 1 
q A R "c 


5 ! v 


£.p P£_ + ^Pp ) 


1 U X n + C 

KrI a --- 3 » 

6 t 1 - u 2 x°t x°; 

q q 


(4.6) 


The integration over the longitudinal fluctuation propa- 

2 2 

gator gives a x /m contribution. 


d. The magn on dr ag term : This represents the effect of 
fluctuation interaction on the RPA term. There are four 
diagrams corresponding to this term (Fig.23b,c,d and e) . 
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. I ijj 

Y ) 
or d 


e k 


l k a G j<< 4 . G k' + G k«+q»+ G 


k + c] ? f 


U X q> 


i - u x° + 7 

-q 


+ k^q ,i<a G Ji , + 1.' ~5l' * ^ 

2 


(1 - u 2 X* ; X~ , ) ( 1 - U X°^7) 


We substitute these into eqn. (2.1). l\s usual consider 
C['=0 in the electron propagators. We will then have terms 
of the form 


~£ q a k a ( + A 3 i^' + AC k f) x (a k + AB k + A C ld q fi k 


k' ''6 


1 I 

P 3_' 


U 3 X 


q’ 


( )( ) 

= l' (q.kM f k , l (q.k) f k l r L (-qJ) f 1 ( q ' ) 


We notice that on K or !<' 

I ( q . t) % (scalar) = D. 
k k 

“ ^ G. , . G. . G. . . 

3 v) j< T j< f is t 
V 1 


summation it will vanish, 
Here we h nv e defined 


“ ^ a r'i + APj r<’ + 


A 2 C^, ) 


si nc e 


So nothing interesting survives if we j ut 3 _ r -w. Now let 
us consider the full expression. This will involve the 


factors of the form 
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B G k' + / + G k'+a‘+ 


To the lowest order in ( q . q ? ) the term can be written as 


A j + A t> 


L 1 ) ( q. q ' ) -t- 0( | n , q 1 { ' ) 


We procaed as before and find (after collecting a, b, c, 
d &, e „ terms and analyzing) that the surviving term is 

(mDq 2 )" ,D S’ . iSl. 1 4 2 ( 4 I ■ ” T - 1 - 


o - q 1 c\ r , r , ) 
P ‘a '3 C|' -q’ 


(4. 7) 


The most dominant term will be that from the pole of T 


T 


and will contribute a m(T) T 


5/2 


term at low temperatures 


and hence negligible companion to T etc. terms. 


Collecting all the thermal contributions we can write 
the expression for D(T) as 

-2 


n(T , .„ T , r -1 n RPA , t- x . 3/2 

D(T) = l'l(T) [ii -0 - A ~~7 ~ 3 in (w 

m 

-Cm (t/m) ^ ] • 


(4. 8) 


The coefficients A, B and C are given in terms of single 

particle density of states and its Derivatives at fermi 

energy. hot ice that at low temperatures there is a 

2 

strong T dependence, while at high temperatures D(T) 
is expected to decrease as xm(T). 



5. 


D is cu ss io n 


The main qualitative features of our results Eq. (4.8) 

, 2 

for D are: 1 ) an enhanced T dependence (enhanced doubly in 

- 1 3/2 

powers of rti , a large quantity) , 2) h reduced T depen- 

3 /2 

dence (reduced by one power of m) . The presence of a T 

conflicts with the phenomenological result of Izuyama and 

5/2 

Kubo who find only a T term. It also conflicts with 

2 8 

the result of Corrias and de Pasquale who calculate 

(for a strong ferromagnet <n^> = w , <n^> = 0) the spin 

wave spin wave vertex C(k ? k ! ) (see Eq. (1*9)) using a Ward 

2 2 

like identity and show it to be of the form k k } for small 

2 

k ? k ! . (Both those authors find a T term.) We have no 

3 /2 

good explanation for this discrepancy. The T term 
arises from a self energy correction to the single particle 
propagator in the particle hole term . We have 

looked for and not found any cancellation of this term, 
to leading order in T. We are unable to cast our result 
in the I zuy am Kubo form of i.e., we do not write the total 
energy of the system as a functional of the number of spin 


w aves 


and of electrons. It is not clear that this can be 


do ne . 



CHAPTER VII 


CONCLUSION 


In this chapter we make a few concluding remarks. 

We briefly summarize cur approach and results, point out a 
few areas where more work needs to be done, and briefly 
discuss other phase transitions in fermi systems where 
fluctuation effects of the type considered might be relevant. 

In this thesis, nearly and weakly ferromagnetic 
fermi systems have been considered. It is found that the 
thermal or temperature dependent properties are dominated 
by spin fluctuations and their interaction. Spin fluctua- 
tion excitations are of low energy, and are fairly stiff 
(their energy increases only slowly with increasing 
fluctuation wavevector q) . They couple strongly amongst 

themselves and the coupling is basically of short range 
-1 

('Vkp ) # The characteristic temperature dependences 
arising in physical properties as a result of coupling 
between these spin fluctuations can be explicitly deter- 
mined in two extreme physical regimes i.a. the quantum 
or low temperature regime and the classical regime. 
Expressions retaining the leading temperature dependent 



terms can be written down which cover these extremes and 
the intermediate region as well. Because the systems is 
by assumption at temperatures far below the free fermi 
temperature T ° 9 effects due to more than two correlated 
spin fluctuations are shown to be small. Since the two 
correlated spin fluctuation term behaves very much like 
the one spin fluctuation term* a simple Hartree like 
approximation (self consistent mean fluctuation field 
theory) is adequate generally. These c o ns id er at io ns 
break down near the critical point. This critical regime 
shrinks exponentially as T^-^0. Because of the inter- 
mediate coupling nature of the problem, while the tempera- 
ture dependence can be found, its magnitude (i.e. the 
coefficient) cannot be determined accurately from theory. 

We show that these intermediate coupling effects determine 
weakly temperature dependent or zero temperature quantities 
like the fluctuation coupling vertex or zero temperature 
susceptibility. These are either taken from experiment or 
parametrized in a simple way. The results obtained from 
this approach are presented in Chapters II to VI, In 
Chapters II and III, we compare our results with experi- 
ment for x(T) and C (T) of liquid He 3 . The fluctuation 
coupling vertex is parametrized by a srmple step function. 

We now indicate some areas of incompleteness. In 
the nearly ferromagnetic regime, the specific heat has been 



calculated in the quas Lharmonic approximation (QHA). We 
have not investigated in detail the error incurred in 
making this approximation, i.e. the size and temperature 
dependence of the leading correction term. Such an estimate 
was made for X( T ) in Chapter II, which led to the conclusion that 
the fluctuation correlation corrections are small for 
(T/T°) — x < < 1 . In the case of free energy, some simple 

estimates show thax the diagrams left out in UHA‘ do not 
contribute to leading order. However, this point needs 
more careful scrutiny. From experience with systems under- 
going phase transitions, one finds that ^ HA (which trans- 
lates into mean field theory) is quite good away from the 
critical region. Since here one is almost but never in 
the critical region, we expect cur results to be realistic. 

In the paramagnetic (T > T ) regime, an important 
feature is the observation of spin wave like modes. 

These seem to be well defined for intermediate values of 
q, to be almost as ? stiff 1 as in the ferromagnetic phase 
for these values of q, and to become overdamped beyond a 
fairly sharply defined cutoff q . To understand this 
phenomenon in terms of spin fluctuations and their inter- 
actions; one has to investigate in detail the propagation 
of a spin triplet electron hole pair in the thermal, 

classical spin fluctuation field. This has been done 

79 80 _ . 81 
recently by Ihoriya, Klenin and Hertz, and R amakr is nnan. 



We do not discuss this effect here; its influence on 
results obtained here for nearly ferromagnetic systems is 
expected to be small. 


In the ferromagnetic regime, we have not made any 
detailed comparison between theory and experiment, either 
for magnetization m(T) or for spin wave stiffness D(T). 

The main problem is the band structure of the weak ferro- 
magnets. Their electronic energy spectrum is not known 
and is perhaps not well approximated by a free electron 
like form. However, m(x)/m(0) and D(x)/D(0) curves, plotted 
as a function of x for various values of x , would be of 
value. Obtaining this requires one to solve simultaneously 
and self consistently for m(x) and the fluctuation ampli- 


tudes D (£ ), D (5 ). One expects a nearly ’universal' 

o o 0 9 

m (x) 

curve because while at low temperatures = 1 - A (~j) , 

m 2 (t) T < 

at high temperatures (near x ) — = (1 - B — ) . 

m (0) c 

A and B depend on the fermi energy. In the case of Arrott 

2 

plot of m (x) vs (H/m), the deviations from linearity are 

2 , 

directly expressed in terms of m lx). The deviations 
observed are of the right sign and show the right trend 
with t j a more detailed comparison is in order. In case 
of D(x), there are no data available for weak f e rr om agnet s . 


We have no satisfactory understanding of the 


discrepancy between the T 


3/2 


term present in our theory 
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5/2 

and the T term obtained by others for D(T). An alter- 
native way of calculating D(T) is being attempted. This 
is basically a time dependent Ginzburg Landau (TDGL) theory, 

and has been used recently for example, in discussion of 

82 

dynamic critical phenomena. At low temperatures there are 

quantum effects, i.e. fluctuations with Matsubara 

qrn 

frequencies m/0 need to be considered also. Thus the usual 
TDGL theory, which deals with classical order parameter 
fluctuations, needs to be generalized. In this semipheno- 
menological formulation one deals directly with the 
equation of motion of the magnetization l v i in terms of 
functions of M and V |V|. The formal machinery of Chapter VI 
is circumvented, and the formulation is physically more 
d ir ec t . 

Are the kind of fluctuation effects considered here 
important for ferrni systems exhibiting other phase transi- 
tions? In the superconducting phase transition, because 
of the large zero temperature coherence length £ , and 
because fluctuation coupling is sharply peaked in wave- 
vector and frequency space, ^ fluctuation effects are 
small. The BCS theory which is a generalized H artree-Foc k 
theory similar to that of Stoner and Wohlfarth for ferro- 
magnetism, works quite well. In spin density wave or 

itinerant electron antiferromagnetic phase changes, the 

83 

situation is not clear. Moriya and have shown 



that one expects considerable spin fluctuation effects. 

The experimental evidence can be interpreted variously, 
and strong band structure effects are passably present 
in , s ay, C r . 

It is possible that the recently observed CDW phase 

transitions in layered c h alco ge nid es ^ 9 ^ may afford another 

example where fluctuation effects of the type considered 

here contribute significantly to temperature dependent 

properties. One has, here too, an electron hole in a 

fermion system with strong short range coupling. Experi- 

2 

mentally, for example, a large I term in the resistivity 
is seen over a wide temperature range. The size is far 
too large to be explained by a simple Baker scattering 
mechanism. On the other hand, if the resistivity is due 
to scattering from density fluctuations with q ^ 

is proportional to their number, which, in the low tempera- 
ture range, is doubly enhanced (Chapters II and V), i.e., 

pf 1 uc t . ^ (t /T ) a The resistivity also levels off to a 
nearly linear form at higher T 0 



APPENDIX I 


In this Appendix we write the Hamiltonian of an 
interacting Fermi system with short range forces (zero 
range for simplicity) in two rotationally invariant forms. 
The interaction term is in the form of interacting spin 
and density fluctuations. Though the forms look different, 
a proper resummation of diagrams or perturbation series 
leads for example to the same RPA like criterion for ferro- 
rnagnet ism . 


Consider the Hubbard Hamiltonian 


H = KE 4 - ) Un. n. 

1 if x 4/ 


( 1.1 ) 


the Stoner Hamiltonian being of course just a continuum 
form of 1.1. We shall work with 1.1; all the results 
obtained here are valid for the Stoner Hamiltonian too. 

The second term in 1.1 just represents the repulsion when 
two fermions are on the same site. Because of the exclusion 
principle, they are necessarily of opposite spin. We can 
write 


Un... n, , = ~ U ( n ^ + n 4 J -jS..S 


9 ^ 

_ f S ..5 . 

if iv 3 i-i 


' if ‘ ' i V 2 

- T U( n, +n . - 


( 1. 2a) 


1 \ 2 


4 


if' xV 


( 1 .2b) 
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In ( 1 .2 a) , except for a one body term, which can be absorbed 
by redifining the K.E. term, we have a quadratic form involv- 
ing spin fluctuations, i.e. 


U n . n . . 
r't- xit 


j- U ( n . +n . ) 
1 r i k 


I 4s 



(1.3) 


This rotationally invariant quadratic form is quite suitable 
for the functional integral transformation. .This transforma- 
tion is described in Chapter II, the diagrammatic rules being 
given in Appendix II. 


It may appear from these rules and the resulting 
diagrams that the criterion for divergence of X n i 3 1-U’fg =C 

CJ U p 

where U'=~. (Ws calculate the spin fluctuation propagator - 

with the diagram — • ^ for the self energy. 

Here straight lines — ^ — are electron lines and ~^\ ’*’ S 

the bare vertex. The bare fluctuation propagator is denoted 

by a dotted line .) Thus the criterion appears to differ 

from the RFA. Actually one can obtain the RPA in our dia- 
grammatic formalism with the Hamiltonian (1.3) functionally 
transformed. We need to sum the series 



for the fluctuation self energy. This leads to the criterion 



= 0 . 


(1 


0 


i.e. (1 - 


l 1.4) 



We assume that this is done, where necessary, in our work 
in the body of the thesis. 

The form (1. 2b) of Lin. ,n. , is ani'chor rotationally 

1 \ 1 vV 

invariant quadratic form. In this both density and spin 
fluctuations appeir qu adr at ic ally . This form .is a general 
two body interaction, involving both spin end density fluctua- 
tions. We use this form in Chapter III. it can be used else- 
where also. Again we remark that the RPA or Stoner criterion 
for ferromagnetism can be obtained (after the functional 
transformation) by a proper resummation of diagrams for 
fluctuation self energy. 

The variety of forms discussed above is already known 

9 1 

for the one impurity problem (Anderson model) . Keiter has 
shown there how resummation of diagrammatic perturbation 
theory terms removes apparent differences. More generally, 
the differences are an expression of the tact that we have 
an intermediate coupling problem where interaction vertices 
cannot be determined accurately. Retaining another set of 
diagrams would lead to a different U eff in, say the Stoner 


c r iter io n. 



APPEWDIX II 


We give h _-re the rules for drawing various diagrams 
involving interaction between the fer.nin.i and the fluctua- 
tion fields. The Fourier transform of trie spin density 
operator is defined as 
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The transverse fluctuation creation operators are defined 
as 
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IMow from the equation (2 a. 9) of Chapter II, the 


vertices are given by 
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A factor of (-1) for each closed loop and a factor 
1/n! for the order n of the diagram are introduced. The 
combinatorial factor, after considering ail topologically 
identical diagrams, can be written as (1/5) where 5 is 
the symmetry numb.r of the diagram (T. ii.fi ice, Ref. 45). 

A summation over all internal variables is to be performed. 


For the magnetic phase of the system, the vertex for 
coupling to the molecular field is 


/ * 



O 2(^) 1/z 


1 + 

2 a k + ua a k d* 


The effect of the molecular field is cc change the single 

1/2 -f 

particle energy to E^ = + (U/j3) ? 0 and the 

O 1 / 9 

propagator G^ to G^ a = ( + (U/g) ‘~% q ). 
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